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ABSTRACT 


The parameters of Low Probability of Intercept (LPI) radar signals are hard to 
identify by using traditional periodogram signal processing techniques. Using the 
Wigner Distribution (WD), this thesis examines eight types of LPI radar signals. Signal 
to noise ratios of 0 dB and -6dB are also investigated. 

The eight types LPI radar signals examined include Frequency Modulation 
Continuous Wave (FMCW), Frank code, PI code, P2 code, P3 code, P4 code, COSTAS 
frequency hopping and Phase Shift Keying/Frequency Shift Keying (PSK/FSK) signals. 
Binary Phase Shift Keying (BPSK) signals although not used in modem LPI radars are 
also examined to further illustrate the principal characteristics of the WD. 
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I. INTRODUCTION 


A, LPI RADARS 

On the modem battlefield, radar systems faee serious threats from EA (Eleetronie 
Attaek) and ARMs (Anti Radiation Missiles). To perform the neeessary target deteetion 
and traeking and simultaneously hide themselves from an enemy attaek, radar systems 
should be Eow Probability of Intereept (EPI). [Ref. 1] 

Eow probability of intereept radars have the property of low power, wide 
bandwidth, frequeney variability or other design features to make them diffieult to be 
deteeted or identified by passive intereept reeeiver deviees sueh as Eleetronie Support 
(ES) or Eleetronie Intelligenee (EEINT) reeeivers. [Ref. 2] EPI radars attempt to provide 
deteetion of targets at longer range than intereept reeeivers ean aecomplish deteetion of 
the radar. The sueeess of EPI radars is measured by how hard it is for them to be 
deteeted. 

The EPI radars typieally have the following features: 

(1) Very low peak power (on the order of IW) 

(2) Radiated energy spread over a wide angular region, over a long time interval, 
and over a wide frequeney band. 

(3) Continuous Wave (CW) radiation with a large time bandwidth produet. 

(4) Eow sidelobe transmit antenna and earefully seleeted earrier frequeney. 

(5) Redueed reeeiver noise temperature and loss. [Ref 1, 3, 4] 
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The transmitter makes use of sophisticated frequency and phase modulation to 
spread the signal bandwidth making the signal hard to intercept. The receiver makes use 
of the appropriate matched filter so that the radar performance is similar to that of 
traditional pulsed radar radiating the same amount of average power. 

In this thesis, Binary Phase Shift Keying (BPSK), Frequency Modulation 
Continuous Wave (FMCW), polyphase code modulation (Frank, PI, P2, P3, P4), 
COSTAS frequency hopping, and PSK/FSK signals will be analyzed using the Wigner 
distribution in order to examine detection capability and to extract the modulation 
parameters. 

B, TIME-FREQUENCY ANALYSIS (TEA) DETECTION TECHNIQUES 

In many operation radar environments, the non-stationary nature of the received 
radar signal mandates the use of some form of time-frequency analysis (TFA), which can 
clearly bring out the non-stationary behavior of the signal. The important virtue of TFA 
is that it provides an indication of the specific times during which certain spectral 
components of the signal are observed. [Ref. 5] 

Time-frequency analysis allows time and frequency resolution to be controlled 
independently, but when used for analysis of multi-frequency component signals this 
approach is vulnerable to cross-terms arising midway between frequency components. 
There are five techniques, which are in current use. [Ref. 6] 

1. Short-Time Fourier Transform (STFT) for a discrete data x{l) is defined as 

00 _2 -t ^ 

STFT(l,k)= ^ (/-«)• " (1.1) 
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Estimating of the TEA using the STET involves computing a Eourier transform (ET) of a 
discrete time series x{l) using a sliding window . At each time index /, a A point ET 
is computed, giving a spectral estimate at the frequency index k. The window length N 
is chosen to optimize the time and frequency resolution; the signal is assumed to be 
stationary within the window. Signal energy density can be estimated from the squared 
modulus of the STET. This is called the spectrogram. [Ref. 6] 

SP{l,k) = \STFT{l,kf (1.2) 

2. Wigner Distribution of a continuous signal x{t) is defined as 

00 

WD(t,Q))= f x{t + -)-x\t--)-e^^^dT (1.3) 

•'2 2 

r=—00 

where co is angular frequency In f , the integral is from -go to oo , and the * indicates the 
complex conjugate of the signal x{t) . Various formulations of the WD are possible for 

sampled time series. One of these uses a moving window with 2A-1 nonzero values. 

[Ref 6] 

3 . Wigner Distribution (WD) for a discrete data x{l) is defined as 

WD(l,k) = 2 Wj^(n) ■ x(l + n) ■ X* (I - n) ■ e ^ (1-4) 

n=-N 

The WD offers very good time and frequency resolution compared with the SP, but the 
cross-terms generated from analysis of multi-component signals can make visual 
interpretation of the TEA difficult. If another window function with 2M -1 nonzero 
values is introduced, the cross-terms can be reduced. [Ref 6] This gives the smoothed 
WD. 
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4. Smoothed Wigner Distribution (SWD) for a discrete data x{l) is defined as: 

SWD{l,k) = 2 w^{n)\ ^ •x(/ + m + n)-x (/ + m-«) >-e ^ (1-5) 

n=-N \^m=~M J 

The discrete time of WD (1.4) is periodic in n , but for SWD the discrete time is periodic 
in In . So unless the signal x{l) has been sampled at least twice the Nyquist rate, there 
will be aliasing in the WD. This problem can be solved either by sampling the signal at 
twice the Nyquist rate or by estimating the TFA of the real valued time series x(/) . [Ref 
6 ] 

5. Choi-Williams Distribution (CWD) for a discrete data x(/) is defined as: 

N-\ I M-1 

CWD{l,k) = 2Y,w,{n)\ X 

n=-N m=-M 

( 1 . 6 ) 

where the parameter cr can be adjusted to attenuate cross-terms. The CWD is a further 
refinement of the SWD, in which extra smoothing is achieved with an exponential 
weighting function. Lower values of cr attenuate cross-terms but degrade time and 
frequency resolution. For large values of cr, cross-terms are not suppressed at all, and 
the CWD becomes equivalent to the SWD. A good choice of cr is in the range 0.1 to 10. 
[Ref 6] 



In the following chapters, the WD and its results will be introduced and examined. 
The SWD and the CWD are for the further studies of the TFA techniques. They will not 
be examined in this thesis. 
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c. 


PRINCIPAL CONTRIBUTIONS 


The objective of this thesis was to develop the Wigner distribution (in 
MATLAB® [Ref 7]) and use it to analyze and extract the parameters from a variety of 

complex LPI radar modulations. In order to use the Fast Fourier Transform (FFT), a 
special kernel transformation is developed and two simple examples are worked in order 
to insure that the kernel transformation is calculated correctly. 

The 114 signals analyzed include BPSK, FMCW, Frank code, PI code, P2 code 
P3 code, P4 code, COSTAS, PSK/FSK with frequency hopping and two test signals. 
Using the Wigner distribution, this thesis successfully extracted all the parameters of the 
LPI signals except the PSK/FSK with target frequency hopping. 
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D, THESIS OUTLINE 

Chapter II presents a complete introduction of the Wigner Distribution (WD). 
Two simple examples are provided in this chapter. 

Chapter III presents the Binary Phase Shift Keying (BPSK) signals and the 
analysis of BPSK signals using the Wigner distribution. 

Chapter IV presents the FMCW signals and the analysis of FMCW signals using 
the Wigner distribution. 

Chapter V presents the polyphase signals (Frank code, PI, P2, P3 and P4) and the 
analysis using the Wigner distribution. This chapter also presents the comparison of 
polyphase signals. 

Chapter VI presents the COSTAS frequency hopping signals and the analysis of 
COSTAS signals using the Wigner distribution. 

Chapter VII presents the Phase Shift Keying/Frequency Shift Keying (PSK/FSK) 
combined with frequency hopping signals and the analysis of PSK/FSK signals using the 
Wigner distribution. 

Chapter VIII presents some concluding remarks and recommendations. 

The Appendix contains the MATLAB® m-fdes of the Wigner distribution and the 
LPI radar signal matrix. 
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II. THE WIGNER DISTRIBUTION 


The Wigner Distribution (WD) introduced by Wigner in 1932 as a phase 
representation in quantum mechanics gives a simultaneous representation of a signal in 
space and spatial-frequency variables. The Wigner distribution has been noted as one of 
the more useful time-frequency analysis (TFA) techniques for signal processing. 

A, WIGNER DISTRIBUTION EQUATION 

Recall (1.3), the Wigner distribution of a continuous input signal x{t) is defined 
as 

00 

W{t,co) = I x{t + Tl2)-x*{t-Tl2)e^^‘"''dT (2.1) 

-00 

where t is the time variable and o) is the frequency variable. The Wigner distribution is 
a two-dimension function describing the frequency content of a signal as a function of 
time. 

This continuous time and frequency representation can be modified for the 
discrete sequence x{l), where / is a discrete time index, / = ...,-1,0, 1, .. .The discrete 
Wigner distribution (WD) is defined as 

00 

W{l,a)) = 2Y,x{l + n)x\l-n)e-^^‘^" ( 2 . 2 ) 

«=-oo 

Further modification results in the pseudo-Wigner distribution or windowed-Wigner 
distribution, which is defined in (2.3), the same as (1.4). 

N-l 

= 2 ^ x{l + n)x*{I - n)w{n)w{-n)e^^^""' (2.3) 

n=-N 


1 



where x(/) is a diserete input signal with / from -go to oo , w{n) is a length 2N real 
window funetion with w(0) = 1. Here N must be as large as possible within the limits of 
an aeeeptable eomputational eost beeause a large N gives more output samples, yielding a 
smoother result. Considering the reetangular window funetion with a magnitude equal to 
one the WD beeomes [Ref 8] 

V-l 

r(/,®) = 2 ^ x{l + n)x{l-n)e-^^‘"" (2.4) 

n=-N 

Using fi{n) to represent the kernel funetion 

f(n) = x(l + n)x*(l-n) (2.5) 

the WD beeomes 


N-l 

W{l,co) = 2j^ 

n=-N 


where the eontinuous frequeney variable o) is sampled by [Ref 9, 10] 


( 2 . 6 ) 



(2.7) 


and where A: = 0,1, 2,...2N- 1. The kernel indexes are modified to fit the standard 
Diserete Fourier Transform (DFT). Sinee 

f{n) = fX-n) (2.8) 


the kernel is a symmetrie function, so the Discrete Fourier Transform (DFT) of the kernel 
is always real. 

From equation (2.6) and (2.7), the WD becomes 
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(2.9) 


r(/,^) = 2y/,(„)exp(-.^) 
2N 2N 


Adjusting the limits of n in order to use the standard FFT algorithms, (2.9) beeomes 


»r(/,^) = 2y/(„)exp(-.^) 
2N ^ 2N 


( 2 . 10 ) 


Note that in (2.10) the kernel funetion has been adjusted to f, (n) , where 


// («) = 1 



[Mn-m 


0< n< N 
n = N 


N + \<n<2N-\ 


( 2 . 11 ) 


The resulting WD is 


27^-1 

W(l,co) = 2Y,fi(n)e-^^-" 

n=0 


( 2 . 12 ) 


2^1 , 

W{l,k) = 2Y, f{n) exp(-7^) 

N 


(2.13) 


Equation (2.13) is the final WD equation, and will be used on the following paragraph to 
show the eomputational proeedure of the WD values. 


B, TWO SIMPLE EXAMPLES 

1, Example With Real Inputs 

Suppose there is an input signal 


x(/) = 12,4,3,6,1,7} 

= {2(/ = -3), 4(1 = -2), 3(1 = -1), 6(/ = 0), 1(1 = 1), 7(1 = 2)} 


(2.14) 
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where the length of the input signal x(/) is 2N= 6,oxN=3. / is a diserete time index 


from -N to N-\. Note that x = 0 for / < -4 or / > 3 . Reeall (2.11), sinee now A^=3, so 



f^{n) , 0<n<2 
0 , n = 3 


4<n<5 


(2.15) 


For example, for / = -3, the kernel f_^ {n) beeomes 

U{n) = {U{n = 0),/_3(n = l),/_ 3 (n = = 2),f_,{n = A),U{n = 5)}, from (2.15) 


/_ 3 '(n = 3) = /_ 3 (n = 3) = 0, f_;{n = A) = U{n = A-6) = f_,{n = -2) and 


/_ 3 '(n = 5) = /_ 3 (n = 5-6) = /_ 3 (n = -l) 

Henee, 

U{n) = {/_ 3 (n = 0 ),/_ 3 (n = l),/_ 3 (n = 2)AU{n = -2),f_,in = -1)} 


Reeall that f,{n) = x{l + n)-x\l-n), 


and x(/) = {2(1 = -3), 4(1 = -2), 3(1 = - 1), 6(1 = 0), 1(/ = 1), 7(1 = 2)} 

and x(l) = 0 for / < -4 or / > 3 . 

f_j(n) for input signal x(/)is eomputed as follows 


/_3 (n = 0) = x(-3 + 0) • x’ (-3 - 0) = x(-3) • x* (-3) = 2-2 = 4 
/_3 (n = \) = x(-3 +1) • X* (-3 -1) = x(-2) • x* (-4) = 4-0 = 0 
/_3 (n = 2) = x(-3 + 2) - X* (-3 - 2) = x(-1) - x* (-5) = 3-0 = 0 
/_ 3 (n = 3) = 0 

/_3 (n = 4) = /_3 (n = -2) = x(-3 - 2) - x* (-3 + 2) = x(-5) - x* (-1) = 0 - 3 = 0 
/_3 (n = 5) = /_3 (n = -\) = x(-3 -1) - x* (-3 +1) = x(-4) - x* (-2) = 0-4 = 0 
So, /_;(«) = 14,0,0,0,0,0} 
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And similarly for / = 0, 


/g (« = 0) = x(0 + 0) • X* (0 - 0) = x(0) • x*(0) = 6 ■ 6 = 36 
/g {n = \) = x(0 +1) • /(O - 1) = x{\) • /(-I) = 1 ■ 3 = 3 
/o (n = 2) = x(0 + 2) • X* (0 - 2) = x{2) ■ x (-2) = 7-4 = 28 
/o(n = 3) = 0 

/o(^ = 4) = /o(« = -2) = x(0-2)-/(0 + 2) = x(-2)-/(2) = 4-7 = 28 
f^{n = 5) = /o(n = -1) = x(0 - 1) • x’(0 + 1) = x(-l) • / (I) = 3 • 1 = 3 
So, /o'(n) = {36,3,28,0,28,3} 

Repeating the above procedures, the kernel matrix for / = -4 to 3, and n = 0 to 5 is as 
shown in Figure 1. 


/=3 

0 

0 

/ 


0 


0 

\ 

1=2 

49 

0 


0 

^minet 

ric^ 


/=l 

1 

42 

0 ^ 

► 


0 

\ 

42 

/=0 

36 

3 

28 

0 

28 

3 

/=-! 

9 

24 

2 

0 

2 

24 

/=-2 

16 

6 

0 

0 

0 

6 

/=-3 

4 

0 

0 

0 

0 

0 

/=-4 

0 

0 

0 

0 

0 

0 


« =0 

«=1 

n=2 

n=3 

rt=4 

n=5 


Figure 1. The Kernel fi {n) Matrix for the Real Six Input Example. 
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The second step is to use (2.13) to calculate the Wigner distribution. As an example of 
the calculation, one can pick any / and k to examine the values inside the WD matrix. 

For example, choose / = 1, A: = 2 and A= 3. The WD is 

. Trkn 

W(l = l,k = 2) = 2Yf, (n)exp(-y—) 

N 

2 - 3-1 , 2n 

= ) (2.16) 

n=0 ^ 

^ 1 ^ 77' T1 

n=0 2 ) 

From the kernel matrix in Figure 1, the kernel function for / = 1 is 
(n) = {1,42,0,0,0,42}. Which means that 
/'(O) = 1,/;(1) = 42,/'(2) = 0,/'(3) = 0,/'(4) = 0,/'(5) = 42 . 

So from (2.16), the WD for / = 1, A: = 2 is 

ir(l,2) = 2y(0).exp(-y^) + 2/;'(l).exp(-y^) + 2/;'(2).exp(-y^) 

.2y;(3).exp,-y^).2y;,4).exp(-y^).2/,5).exp,-y^) 

= 2•l.(0) + 2.42.(-0.5000-0.86600 + 2-0 + 2-0 + 2•0 + 2.42.(-0.5000 + 0.8660i) 
= -82 

Repeating the above procedures, gives the WD matrix at discrete time index / = -4 to 3 
and discrete frequency index A: = 0 to 5 as a symmetric matrix shown in Figure 2. An 
important feature in this Wigner distribution computation result is that all the components 
in WD matrix are real. 
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/=3 

0 

0 

0 

0 

0 

0 

/=2 

98 

98 


98 

mmeti 

9 

\D 

00 

^98 

/=1 

170 

8^ 

j 

00 


X-82 

▼ 

/=0 

196 

22 

10 

172 

10 

22 

/=-l 

122 

62 

-34 

-70 

-34 

62 

/=-2 

56 

44 

20 

8 

20 

44 

/=-3 

8 

8 

8 

8 

8 

8 

/=-4 

0 

0 

0 

0 

0 

0 


k=0 

k=\ 

k=2 

A=3 

k=A 

k=5 


Figure 2. The WD Matrix W(l,k) for the Real Six Input Example. 

2, Example With Complex Inputs 

Consider the input signal 

x = l + j-Q 


where 

/ = cos(2;r fj) 
Q = sin(2;rX0 


(2.17) 


(2.18) 

(2.19) 


If carrier frequency /,=lkHz; sampling frequency f^ = 7kHz; and t = 0: 


—: —, then the 

fs fs 


first eight input points from discrete time index / = -4 to 3 is 
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x(/) = {l + Oz, 0.6235 + 0.7818/, -0.2225 + 0.9749/, -0.9010 + 0.4339/, 
-0.9010 - 0.4339/, -0.2225 - 0.9749/, 0.6235 - 0.7818/, 1 + 0/} 


Consider the value when / = 0, n = 3. 

Using (2.11) with an input length 2N = 8 or = 4. The kernel is 





,0 < « < 3 
,n = 4 
,5<n<l 


( 2 . 21 ) 


Here // {n) = {f, (1), / (2), f, (3), 0, / (-3), f, (-2), /(-!)} 

Sinee f, (n) = x{l + n)-x*{l-n) , the kernel at / = 0, n = 3 is 

/o' (3) = /o(3) = x(0 + 3) • /(O - 3) = x(3) • x*(-3) 

= 1-(0.6235+ 0.7818-/)* 

= 0.6235-0.7818-/ 

Repeating the same proeedures as discussed in real input case, the kernel matrix for the 
eight inputs complex input example is shown in Figure 3. 
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3 


1 

0 

0 

0 

0 

0 

0 

0 

1 

- 0.2225 + 
0.97491 

0 

0 

0 

0 

0 

- 0.2225 - 
0.97491 

0.9999 

- 0.2226 + 
0.97491 

- 0 . 9010 - 

0.43391 

0 

= 0 ^ n=l 

0 

1 

i 

0 

- 0.9010 + 
0.43391 

- 0.2226 - 
0.97491 

1.0001 

- 0.2225 + 

0.97491 

- 0.9009 - . 
0.43391 ( 

^^ 623 ^ 

0.78181 

\ 

0.6235 + 

0.78181 

- 0.9009 + 

0.43391 

- 0.2225 - 

0.97491 

1.0001 

- 0.2225 + 
0.97491 

- 0.9009 - 
0.43391 

0.6235 - 
0.78181 

0 

0.6235 + 
0.78181 

- 0.9009 + 
0.43391 

- 0.2225 - 
0.97491 

0.9999 

- 0.2226 + 
0.97491 

- 0 . 9010 - 

0.43391 

0 

0 

0 

- 0.9010 + 
0.43391 

- 0.2226 - 
0.97491 

1 

- 0.2225 + 
0.97491 

0 

0 

0 

0 

0 

- 0.2225 - 
0.97491 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

2 

^ n ^ 

5 

6 

7 


2 


1 


0 


/ 


-2 


-3 


4 


Figure 3. The Kernel Matrix for the Complex Eight Input Example. 


Now the WD matrix is examined from Eigure 3, when / = -1, the kernel is 

/_i'(n) = {1.0001, -0.2225 + 0.9749/, -0.9009-0.4339/, 0.6235-0.7818/, 
0,0.6235 + 0.7818/, -0.9009 + 0.4339/, -0.2225-0.9749/} 


Consider the case when I = -1, k= 4. Recall (2.13), the WD for A^= 4 is 

, Trkn 

W{l = -\,k = 4) = 2^ f,{n) exp (-7 -—) 

N 

= (2.23) 

H=0 4 

7 

= 2^/_i'(n)exp(-7n;z-) 

«=0 


Erom (2.22) and (2.23) 
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W{l = -\,k = A) = 2^/_i {n) ■ Qx^^{-jnn) 
«=0 


/i (0) • exp(O) + /i (1) • exp(- 7 ;z-) + fi, (2) • {-jin) + fi^ (3) • 

+/, (4) • + fi, (5) • i-jSn) + fi, (6) • {-jen) + (7) • (-yV;^) 


= 2 - 


1.0001 • (1) + (-0.2225 + 0.97490 ■(-!) + (-0.9009 - 0.43390 ■ (1) 
+(0.6235 - 0.781 Si) • (-1) + 0 + (0.6235 + 0.7818) • (-1) 
+(-0.9009 + 0.43390 ■ (1) + (-0.2225 - 0.97490 ■ (-1) 


= -3.2073 

Again the WD matrix of the eight input points is another real and symmetrie 2Nx2N 
matrix. Note the important feature; The WD is always real whether the input signal is 
real or eomplex. Figure 4 shows the WD matrix of the eight complex input samples. 


2.0000 

2.0000 

2.0000 
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2.0000 

2.0000 

2.0000 
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- 1.3868 

- 2.4943 

2.3923 

9.5036 

7.1225 

- 0.7139 

- 1.8638 

1.7041 

0.3485 

0.0004 

- 1.5882 

12.6307 

6.6749 

- 3.2073 

-^ 

2.1112 

- 1.4231 

0.7965 

0.0004 

- 1.5882 

12.6307 

6.6749 

i- 

( j . 2073 ^ 

=-l,k=4 

2.1112 

- 1.4231 

0.7965 

- 2.4943 

2.3923 

9.5036 

7.1225 

- 0.7139 

- 1.8638 

1.7041 

0.3485 

1.1099 

4.1280 

5.8995 

5.3867 

2.8899 

- 0.1282 

- 1.8997 

- 1.3868 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 

2.0000 


0 i 2 3 4 5 6 7 


3 
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1 

0 

/ 

-1 

-2 

-3 

-4 


Figure 4. The WD Matrix for the Complex Eight Input Example. 
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Figure 5(a) shows the mesh plot of the eight input complex example in time 
domain. The time resolution 1/ is indicated on this plot. Figure 5(b) shows the mesh 
plot in frequency domain. The carrier frequency in this plot is about 900Flz, very close to 
real value IkHz. The frequency resolution 111#samples is also indicated on this plot. 


Time Resolution= l/fs= 1/7000 = 0.1429 ms 

Mesh of8 inputs complex example 



IkHz Carrier frequency Freq. Resolution = fs/2/#samples 



(b) 

Figure 5. WD for Eight Inputs Complex Example (a) 2D Mesh in Time Domain, (b) 2D 

Mesh in Erequency Domain. 
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Figure 6(a) shows the 3D mesh plot of the eight inputs complex example. This 
plot shows the magnitude in both time and frequency domain. Figure 6(b) shows the 
contour plot. The contour plot is a 2D time- frequency domain plot. The magnitude is 
represented by different color as shows in the color bar. 


Mesh ofS inputs complex example 



X 10 '^ Contour ofS Inputs complex example 



Frequency (Hz) i^Hz Carrier frequency 

(b) 

Figure 6. WD for Eight Inputs Complex Example (a) 3D Mesh, (b) Contour. 
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Figure 7(a)(b) shows the 2D time domain and frequency domain mesh plot of the 


101 inputs complex example (Eq. (2.18), (2.19) with t = 0: —;-). The IkHz carrier 

f S -f s 

frequency is more clearly shown in Figure 7(b) than that in Figure 5(b). Comparing with 
Figure 5 and Figure 7 one can find the more input sample points that are used the better is 
the WD output performance. 


Mesh of 101 inputs complex example 



(a) 


Figure 7. 


Mesh of 101 inputs complex example 



Frequency (Hz) 


(b) 

WD for 101 Inputs Complex Example (a) 2D Mesh in Time Domain, (b) 2D 
Mesh in Frequency Domain. 
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Figure 8(a) shows a zoom in of the mesh plot in time domain. The total time 
interval and time resolution 1//^ are indicated on this plot. Figure 8(b) shows the zoom 
in mesh plot in the frequency domain and indicates the frequency interval and frequency 
resolution fJ2/#samples. 


Time resolution = /s 
_ 1 
7000 


Time interval 

= 4^*(#points) 

fs 



Freq. resolution = 2-#s{mples 

= /y/2 = 35 00 Hz 



Figure 8. WD for 101 Inputs Complex Example (a) Zoom In 2D Mesh in Time Domain. 

(b) Zoom In 2D Mesh in Frequency Domain. 
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Figure 9(a) shows the three-dimension mesh plot of the 101 inputs complex 
example. This plot shows the magnitude in both time and frequency domain. Figure 9(b) 
shows the contour plot of the 101 inputs complex example. Contour plot is a 2D time- 
frequency domain plot. The magnitude is represented by different color as shows in the 
color bar. The carrier frequency is also marked on the plot. 


Mesh of 101 inputs complex example 



-50 

0,015 


0.014 

0012 

0,01 

•0 008 

0.006 

0 004 

0.002 


4000 


Time (s) 


0 0 


Frequency (Hz) 


(a) 


Contour of 101 inputs complex example 


I I 



I » 


500 1000 1500 2000 

Frequency (Hz) 


(b) 


10 3000 

IkHz Carrier frequency 


Figure 9. WD for 101 Inputs Complex Example (a) 3D Mesh, (b) Contour. 
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Now consider the WD for two earrier frequeneies with /^j = IkHz , = 2kHz as 


shown in Figure 10. Reeall that / = eos(2;r fj:) and Q = sin(2;r f^t) . Now 


I = cos(2;r f^^t) + oos(2;rg = sin(2;r fj) + sin(2;r 

Figure 10 (a) is the 2D frequency domain mesh plot whieh shows the two earrier 
frequeneies as well as the frequeney eross term. The frequeney of the cross term is 


/cl + fc2 


2 


. Note that the magnitudes' are positive for the two real carrier frequeneies. But 


for the eross term, the magnitude includes both the positive and negative parts. 


Figure 10(b) shows the eontour plot of these two earrier frequencies and their 
eross term. Note that the shape and eolor of the eross term is not like those of the two 
real earrier frequeneies. The eross term eontains many red and blue ‘dots’, but the real 
fc\fc 2 are yellow ‘strips’. 


Figure 11(a) shows the 3D mesh plot. In this plot displays the eross term is 
eombined by many ‘peaks’. The positive peaks represents the red dots, the negative 
peaks represents the blue dots in figure 10(b). 

Figure 11(b) is the 2D mesh in the time domain. It shows the real signal and the 
eross term in the time domain. This plot states that the eross term is eonsisted of a series 
of positive and negative magnitude eomponents in time domain. 


' The label “Magnitude” was used for WD value in all the mesh plots. Note that is 
not a true magnitude. 
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Mesh of two carrier frequencies example 



(a) 


Contour of two carrier frequencies example 



(b) 

Figure 10. WD for Two Carrier Frequeneies Example (a) 2D Mesh in Frequeney 

Domain, (b) Contour. 
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Mesh of two carrier frequencies example 



f, = 2kHz 
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Mesh of two carrier frequencies example 




(b) 


Figure 11. WD for Two Carrier Frequeneies Example (a) 3D Mesh, (b) 2D Mesh in Time 

Domain. 
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C. SUMMARY 

This chapter contains some simple example signals to generate the kernel function 
and Wigner distribution (WD). No matter whether the signals are real or complex, the 
kernel and WD matrix are always real and symmetric. This is a very important feature 
for the Wigner distribution and an important reason why the WD is a good tool to analyze 
signals by the mesh and contour plots. 

In Figures 5 through 11, the carrier frequencies are very easy to identify in these 
simple examples. In the following chapters, the signals contain complex LPI 
modulations and are no longer that simple. The Wigner distribution’s ability to 
determine the phase and frequency modulations for those signals will be demonstrated. 

Figures 10 and 11 illustrate the frequency cross terms and how they are 
represented in the WD mesh and contour plots. It is shown that the cross terms 
sometimes make the signal analysis more difficult. 

In the next chapter, the BPSK signals and its WD analysis will be examined. 
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III. BINARY PHASE SHIFT KEYING (BPSK) 


A, BPSK 

Phase shift keying (PSK) is an angle-modulated, constant-amplitude digital 
modulation. PSK is similar to conventional phase modulation except that with PSK the 
input signal is a binary digital signal and a limited number of output phases are possible. 

Binary phase shift keying (BPSK) has two phases for a single carrier frequency 
(“binary” meaning “2”). One output phase represents logic 1 and the other is logic 0. As 
the input digital signal changes state, the phase of the output carrier shifts between two 
angles that are 180° out of phase. BPSK is a form of suppressed-carrier, square-wave 
modulation of a continuous wave (CW) signal. [Ref 11] 

Figure 12 shows the BPSK modulator phase diagram. Figure 13 shows the BPSK 
modulator output phase vs. time relationship. Normally, there are many carrier cycles per 
modulation period. 


cos (oj: 
(+90°) 

-5 in coj 

(180°) ^ ^ 
Logic 0 

-cos CO J 
(-90°) 


5 in coJ 
^ ( 0 °) 
Logic 1 


Figure 12. BPSK Modulator Phase Diagram. 
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Figure 13. BPSK Modulator Output Phase vs. Time Relationship. 

1. Binary Phase-Coded Pulse Compression 


Phase changes can be used to increase the signal bandwidth of a long pulse for 
purpose of pulse compression. A long pulse of duration T is divided into N subpulses 
each of width 4 . An increase in bandwidth is achieved by changing the phase of each 
subpulse (since the rate of change of phase with time is a frequency). A common form of 
phase change is binary phase coding, in which the phase of each subpulse is selected to 
be either 0 or ;r radians according to some specified criterion. 

The radar receiver correlates the received pulses with the known transmitted 
pulse, usually with a matched filter. Thus output of the matched fdter will be a 
compressed pulse of width and will have peak N times greater than that of the long 

pulse. The pulse compression ratio equals to the number of subpulses N = T/t^ ~ BT , 
where the bandwidth B 
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The matched filter output extends for a time T on either side of the peak response. 
The unwanted but unavoidable portions of the output waveform other than the 
compressed pulse are known as time sidelobes. When the selection of the phase is made 
at random, the expected maximum (power) sidelobe is about UN below the peak of the 
compressed pulse. [Ref 4] 

2. Barker Codes 

One family of binary phase codes that produce compressed waveform with 
constant sidelobe levels equal to unity is the Barker code. Figure 14 illustrates this 
concept for a Barker code of length seven. There are only nine known Barker codes that 
share this unique property. They are listed in Table 1. [Ref 12] 



(b) 

Figure 14. Barker Code of Length 7 (a) A Long Pulse with 7 Equal Subdivisions Whose 
Individual Phase are 0°(+) or 180°(-). (b) Autocorrelation Function for Plot (a). 
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In general, the autocorrelation function (which is an approximation for the 
matched filter output) for code length N, Barker code will be wide. The main lobe 

is 2ti, wide. The peak value is equal to N. There are {N-\)/2 sidelobes on either side of 
the main lobe. Note that in Figure 14(b) the main lobe is 7 and all the other sidelobes are 
one. [Ref 12] 


Code Length 

Code Elements 

Sidelobe Level, dB 

2 

+ -, + + 

-6.0 

3 

+ + - 

-9.5 

4 

+ + -+ ,+ + + - 

-12.0 

5 

+++-+ 

-14.0 

7 

+++—+- 

-16.9 

11 

+ + +-+ — + - 

-20.8 

13 

+++++ — ++_+_+ 

-22.3 


Table 1. Barker Codes. 


The most sidelobe reduction offered by a Barker code is -22.3dB (13 bit code 
length), which may not be sufficient for the desired radar application. However, Barker 
codes can be combined to generate much longer codes. In the following example and 
Figure 15, Barker code B54 is a Barker code combines code length 5 and code length 4. 
But unfortunately, the sidelobes of a combined Barker code autocorrelation function are 
no longer equal to unity. 

B54= { 11101, 11101,00010, 11101 } (3.1) 
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Figure 15. Combined Barker Code with Code Length 5 and 4. 


Figure 16 below shows a basie bloek diagram of the transmitter design. The 
signal x(t) is a CW sinusoid, whieh for this projeet was set to 

x{t) = sin(2;r f^t) (3.2) 

with earrier frequency . The discrete signal xnT is created by sampling x(t) with at a 
sample rate at least twice that of the highest frequency of the CW signal to obey the 
Nyquist sampling theorem and avoid aliasing. The sampling frequency, , is the 

product of the carrier frequency and the sample ratio variable 

SAR=fJl (3.3) 

that is set by the user. To obey the Nyquist sampling theorem, SAR must be greater than 
two to avoid anomalies that may arise from floating point operations right at the Nyquist 
rate. Due to these limitations, the minimum useful sampling frequency is three times the 
signal frequency, or SAR >= 3. [Ref 13] 

Although there are only 13 bits in this Barker code, the actual sequence of bits in 

brkseq also depends on the sampling frequency and the number of complete periods of 
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signal xnT that are to occur within the same timeframe as one Barker code bit. If the 
number of periods per Barker bit (NPBB) is set to one, one full period of the carrier 
frequency fits within one bit of the 13-bit Barker code. 



Figure 16. BPSK Transmitter Block Diagram. 


In this chapter, twelve BPSK signals will be examined. They are numbered from 
one to twelve and listed in Table 2. A IkHz carrier frequency and 7kHz sampling 


frequency are used for these signals. 


No. 

BPSK 

Barker Code 
Length 

NPBB 

SNR 

1 

B 1 7 7 1 s 

7 

1 

Signal Only 

2 

B 1 7 7 1 0 

7 

1 

OdB 

3 

B 1 7 7 1 -6 

7 

1 

-6dB 

4 

B 1 7 11 1 s 

11 

1 

Signal Only 

5 

B 1 7 11 1 0 

11 

1 

OdB 

6 

B 1 7 11 1 -6 

11 

1 

-6dB 

7 

B 1 7 7 5 s 

7 

5 

Signal Only 

8 

B 1 7 7 5 0 

7 

5 

OdB 

9 

B 1 7 7 5 -6 

7 

5 

-6dB 

10 

B 1 7 11 5 s 

11 

5 

Signal Only 

11 

B 1 7 11 5 0 

11 

5 

OdB 

12 

B 1 7 11 5 -6 

11 

5 

-6dB 


Table 2. BPSK Signals. 
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B, WIGNER DISTRIBUTION FOR BPSK 

Consider the same twelve BPSK signals in Table 2 and use them as the inputs for 
Wigner distribution (WD). The WD results of these twelve signals are shown in Figure 
17 through 20. The mesh plots show the frequeney domain of the BPSK signals after the 
WD proeessing. The contour plots show both the frequency domain and time domain of 
the results. 

In Figure 17, the carrier frequency can be clearly identified by the location of 

the highest or lowest peak value in Figure 17(a), or by the center of the symmetric shapes 
in Figure 17(b). 

Secondly, the peak magnitude in Figure 175(a) is about 600. So the 3dB 
bandwidths should be the frequency range between 300 on both sides, which extends 
from 500Hz to 1500 Hz in Figure 17(a). The 3dB bandwidth is also the value of 

divided by NPBB. That is equal to lOOOHz in Figure 17, 18 (NPBB = 1), or 200Hz in 
Figure 19, 20 {NPBB = 5). The 3dB bandwidth is also called the chip rate. 

Thirdly, in Figure 17(a), if we look closely within the 3dB bandwidth, one can 
find there are 15 peaks in the bandwidth range. In other words, there are 14 intervals in 
the range from 500Hz to 1500Hz, and 14 is two times the Barker code length. In Figure 
18(a), 22 intervals are in the 3dB bandwidth range since the Barker code length is eleven. 

Using the same analysis for Figures 17 to 20, we can find the carrier frequency, 
3dB bandwidth and barker code length. However, for NPBB more than one, the analysis 
of Barker code length becomes difficult. 
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Figure 17. WD for BPSK with 7 Bit Barker Code, NPBB = 1, Signal Only, (No.l) (a) 2D 

Mesh in Frequeney Domain, (b) Contour. 
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Figure 18. WD for BPSK with 11 Bit Barker Code, NPBB = 1, Signal Only, (No.4) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 19. WD for BPSK with 7 Bit Barker Code, NPBB = 5, Signal Only, (No.7) (a) 2D 

Mesh in Frequeney Domain, (b) Contour. 
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Figure 20. WD for BPSK with 11 Bit Barker Code, NPBB = 5, Signal Only, (No.10) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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c. 


SUMMARY 


The results of Wigner Distribution analysis for BPSK signals are discussed for 
two situations. When the NPBB (the number of cycles per barker bit) is equal to one, the 
Barker code length can be easily found by the Wigner Distribution. But the disadvantage 
is the carrier frequency information becomes more difficult to determine in low SNR (the 
result figures are not contained in this thesis). Second, when the NPBB is higher than 
one, the carrier frequency and 3 dB bandwidth can be easily found in low SNR, but the 
Barker code length becomes hard to analyze. 

Table 3 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and Barker code bit from Figure 17 to 20. They are done by the 
visible inspection of the WD results. 

In the next chapter, the FMCW signals and their WD analysis will be examined. 
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BPSK 

NO. 

Carrier Frequency 

Bandwidth 

Code Period 

Bits/Code 

1 

100% 

100% 

100% 

100% 

2 

100% 

110% 

100% 

100% 

3 

0% 

0% 

0% 

0% 

4 

100% 

100% 

100% 

100% 

5 

100% 

90% 

100% 

100% 

6 

0% 

0% 

0% 

0% 

7 

100% 

100% 

0% 

0% 

8 

100% 

120% 

0% 

0% 

9 

100% 

150% 

0% 

0% 

10 

100% 

100% 

0% 

0% 

11 

100% 

120% 

0% 

0% 

12 

100% 

150% 

0% 

0% 


BPSKeffectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Table 3. WD Detection Effectiveness for BPSK Signals. 
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IV. FREQUENCY MODULATION CONTINUOUS WAVE (FMCW) 


A, FMCW 

FMCW radar is a technique for obtaining refined range information from radar by 
frequency modulating a continuous signal. The FMCW modulation is readily compatible 
with a wide variety of solid-state transmitters. The frequency measurement of FMCW, 
which is performed to obtain range measurement, can be performed digitally using the 
FFT. The FMCW signals are very difficult to detect with conventional intercept 
receivers. [Ref. 14] 

The frequency modulation used by the radar can take many forms. Linear and 
sinusoidal modulations have been used before. Linear frequency modulation is very 
suitable with the FFT processor to get the range information of targets. The linear FMCW 
emitter uses a continuous 100% duty cycle waveform so that both the target range and 
Doppler information can be measured. The FMCW waveform represents the best use of 
output power available from solid-state devices. Linear FMCW is also easier to 
implement than phase code modulation as long as there is no strict demand on linearity 
specifications over the modulation bandwidth. [Ref. 2] 

The most basic Linear FMCW technique varies frequency in a sawtooth fashion, 
whereby frequency values ramp up over the course of a modulation period and then drop 
back down to the starting frequency to begin the ramp up for the next modulation period. 
Several repetitions of this sequence produce a sawtooth shape in the frequency vs. time 
plot. By processing the reflected waveform, one can extract range information about the 
target. [Ref 15] 
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Figure 21 shows a functional block diagram of a Triangular FMCW transmitter. 
A sinusoidal carrier is frequency-modulated with a message signal that takes the 
triangular shape as shown in Figure 22. The modulated signal is then corrupted with 
channel distortion, which is assumed to be additive white Gaussian noise (AWGN) with 
zero mean and variance. 


Sinusoidal 

Carrier 


Noise 


(AWGN) 



♦ 


FMCW 

Modulator 


>e- 


^ Transmitted 
Waveform 


Triangle 

Generator 


Figure 21. Triangular FMCW Block Diagram. 



Figure 22. FMCW Modulation Signal in Time Domain. 
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The variables in the diagram are as follows. 

= Carrier frequeney 
= Modulation bandwidth 
tm = Modulation period 

The up ramp and down ramp of each triangle must be constructed separately by the 
triangle generator. The modulating waveform for the up ramp can be modeled as 


/m(0=' 


^ AF AF ^ 

F - 1 - 1, 

It 

0, elsewhere 


m 


(4.1) 


Similarly, the modulating waveform for the down ramp is modeled as 


fdit)=- 


. AF AF ^ 
fc~^ j . t, 0<t<t^ 

0, elsewhere 


(4.2) 


In FM systems, the transmit phase (instantaneous frequency) of the carrier is computed as 


^zi(t) = 2;rI / {x)dx 
0 


(4.3) 


Now, assuming the initial conditions of = 0 at t = 0 and substituting(t) and /j(0 
into (4.3) yields the RF phases for the up ramp and down ramp as 


AF AF 1 


2 

AF 


AF 


^^^rf(0 = 2;^-j(/,+—)T-—t k 0 < t < 


(4.4) 


(4.5) 


It is then a simple substitution to arrive at the form of the transmitted signals for 
the up and down ramps [Ref 14, 15] 
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(4.6) 


^.i(0 = ^cos2;r|(X-^)-^ + ^^^|, 
s,,{t) = Acos2n\{f,+^)-t-^t^^, 
^.2 (0 = ^ sin 2;r |(/, - ^) • ^ + ^ ^^ I, 

{ A/t* k'P' 1 

(/c+—)-^- —^ p 


0 < ^ 

0 < ^ 

^ <t <t^ 

0 <t <t^ 


(4.7) 

(4.8) 

(4.9) 


Like those eodes in BPSK, the amplitude A is equal to one. The most important 
three parameters in FMCW code generator are , AF and . 

Twelve FMCW signals will be introduced in this chapter. They are numbered 
from 13 to 24 and listed in Table 4. The “Mod. BW” means the modulation 
bandwidth; 250Hz and 500Hz are used in this table. The “Mod. Period” means the 
modulation period, 20 ms and 30 ms are used here. A IkHz carrier frequency and 
7kHz sampling frequency are used for these signals 


No. 

FMCW 

Mod. BW (Hz) 

Mod. Period (ms) 

SNR 

13 

F 1 7 250 20 s 

250 

20 

Signal Only 

14 

F 1 7 250 20 0 

250 

20 

OdB 

15 

F 1 7 250 20 -6 

250 

20 

-6dB 

16 

F 1 7 250 30 s 

250 

30 

Signal Only 

17 

F 1 7 250 30 0 

250 

30 

OdB 

18 

F 1 7 250 30 -6 

250 

30 

-6dB 

19 

F 1 7 500 20 s 

500 

20 

Signal Only 

20 

F 1 7 500 20 0 

500 

20 

OdB 

21 

F 1 7 500 20 -6 

500 

20 

-6dB 

22 

F 1 7 500 30 s 

500 

30 

Signal Only 

23 

F 1 7 500 30 0 

500 

30 

OdB 

24 

F 1 7 500 30 -6 

500 

30 

-6dB 


Table 4. FMCW Signals. 
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B, WIGNER DISTRIBUTION FOR FMCW 

Consider the same twelve FMCW signals in Table 4 and use them as the inputs 
for the Wigner distribution (WD). The WD results of these twelve signals are shown in 
Figure 23 to 26. The mesh plots show the frequency domain of the FMCW signals after 
the WD. The contour plots show both the frequency domain and time domain of the 
results. 

Take Figure 23 as an example. First, the carrier frequency can be clearly 

found by the location of the highest or lowest peak value in Figure 23(a). In Figure 
23(b), by the location of the arc shape one can also identify the carrier frequency. 

Secondly, the modulation bandwidth AF and modulation period have scaled 
and marked in Figures 23(b) to 26(b) 
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(a) 

contour of 72 ^ 02^5 



Frequency (Hz) 

(b) 

Figure 23. WD for FMCW with AF = 250Hz,t^ = 20ms Signal Only (a) 2D Mesh 

Frequency Domain, (b) Contour. 
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Figure 24. WD for FMCW with AF = 250Hz,t^ = 30ms Signal Only (a) 2D Mesh 

Frequeney Domain, (b) Contour. 
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Figure 25. WD for FMCW with AF = 500Hz,t^ = 20ms Signal Only (a) 2D Mesh 

Frequency Domain, (b) Contour. 
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Figure 26. WD for FMCW with AF = 500//z,t„ = l>0ms Signal Only (a) 2D Mesh 

Frequency Domain, (b) Contour. 
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C. SUMMARY 

For all the LPI signals in this thesis, Wigner Distribution provides better 
parameter estimation for FMCW than the others. Without any other analysis beyond the 
eontour plot, one ean see the most two important parameters, modulation bandwidth and 
modulation period. 

Chapter I introduced the fact that the Wigner distribution is a Time-Frequency 
Analysis (TFA) technique. Since the FMCW signals have the character that both the 
frequency and time were modulated, using the Wigner distribution can yield excellent 
results for FMCW signals, even at SNR equal to -6 dB. See the technical report to be 
published soon for more details. [Ref 16] 

Table 5 below shows the detection effectiveness of the carrier frequency, 
modulation bandwidth and modulation period by visible measurement from Figure 23 to 
26. They are done by the visible inspection of the WD processing results. 

After the BPSK and FMCW signals, the polyphase signals and their WD analysis 
will be examined in the next chapter. 
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FMCW 

No. 

Carrier Frequency 

Mod. BW (Hz) 

Mod. Period (ms) 

13 

100% 

100% 

100% 

14 

100% 

100% 

100% 

15 

100% 

90% 

100% 

16 

100% 

100% 

100% 

17 

100% 

100% 

100% 

18 

100% 

90% 

100% 

19 

100% 

100% 

100% 

20 

100% 

100% 

100% 

21 

100% 

90% 

100% 

22 

100% 

100% 

100% 

23 

100% 

100% 

100% 

24 

100% 

90% 

100% 


FMCW effectiveness on detection 



Parameters 


□ Only signal 
DOdB 

□ (-)6dB 


Table 5. WD Detection Effectiveness for FMCW Signals. 
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V. POLYPHASE 


A, POLYPHASE CODE MODULATION 

Polyphase codes have many useful features such as low range-time sidelobes, 
ease of implementation, compatibility with digital implementation, low cross-correlation 
between codes, the Doppler tolerance of the frequency modulation codes and 
compatibility with bandpass limited receivers. The largest pulse compression ratio that 
can be obtained with the linear Barker codes is only 13 but the Polyphase codes are not 
limited to any finite pulse compression ratio. [Ref 17] 

The phase of the subpulses in Polyphase coded pulse compression need not be 
restricted to the two levels of 0 and n in the binary phase codes. Polyphase codes produce 
lower sidelobe levels than the binary phase codes and are tolerant to Doppler frequency 
shift if the Doppler frequencies are not too large. [Ref. 4] 

The Polyphase code waveforms provide a class of frequency derived phase coded 
waveforms that can be sampled upon reception and processed digitally. [Ref 18] 
Polyphase codes also addresses the reality of the espionage threat as well as the electronic 
attack threat. Programmability and variety are two code characteristics that make 
espionage a more difficult task; thus. Polyphase code will no doubt prove to be the LPI 
waveform of the following decades. [Ref 19] 

In this thesis, the polyphase codes including Frank code, PI code, P2 code, P3 
code, and P4 code are examined. 
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B, 


FRANK CODE 


Frank codes are a family of polyphase eodes that are elosely related to the ehirp 
and have been used suecessfully in LPI radars. The Frank eode has a length , whieh 
is also the eorresponding pulse eompression ratio. If i is the number of the sample in a 
given frequeney and j is the number of the frequeney, the phase of the zth sample ofyth 
frequeney is [Ref 20] 

( 5 . 1 ) 

where i = \:N and j = \:N 


Twelve Frank eode signals are introduced in this ehapter. They are numbered 
from 25 to 36 and listed in Table 6. In this ehapter, the N = 4 and 8 and eyeles per phase 
(CPP) = I and 5 are used. A IkHz earrier frequeney and 7kHz sampling frequency are 
used for these signals. 


No. 

FRANK 

N (Phase) 

Cyeles/Phase 

SNR 

25 

FR 1 7 4 1 s 

4 

1 

Signal Only 

26 

FR 1 7 4 1 0 

4 

1 

OdB 

27 

FR 1 7 4 1 -6 

4 

1 

-6dB 

28 

FR 1 7 4 5 s 

4 

5 

Signal Only 

29 

FR 1 7 4 5 0 

4 

5 

OdB 

30 

FR 1 7 4 5 -6 

4 

5 

-6dB 

31 

FR 1 7 8 1 s 

8 

1 

Signal Only 

32 

FR 1 7 8 1 0 

8 

1 

OdB 

33 

FR 1 7 8 1 -6 

8 

1 

-6dB 

34 

FR 1 7 8 5 s 

8 

5 

Signal Only 

35 

FR 1 7 8 5 0 

8 

5 

OdB 

36 

FR 1 7 8 5 -6 

8 

5 

-6dB 


Table 6. Frank Code Signals. 
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1, Wigner Distribution for Frank Code 

Consider the same twelve Frank eode signals in Table 6 and use them as the 
inputs for Wigner Distribution (WD). The WD results of these twelve signals are shown 
in Figure 27 to 30. The mesh plots show the frequeney domain of the Frank code signals 
after WD. The contour plots show both the frequency domain and time domain of the 
results. 


Take Figure 27 as an example. First, the carrier frequency can be clearly 
found by the location of the highest or lowest peak value in Figure 27(a). Secondly, in 
Figure 27(b), the 3dB bandwidth can be measured as lOOOHz. Since 


h cpp 


(5.2) 


so, the cycles per phase {CPP) can be computed by CPP = 


L 

B 


1000 

1000 


For Frank code 


7 

T = N\=— (5.3) 

since the code period (7) is measured as 0.016 in Figure 27(b), so the phase length (N^) 
can be calculated from = B-T = 1000• 0.016 = 16. The slope which measured in 
Figure 27(b) as 1.6x10 sec^ can be also use to determine the and CPP by the 
following equation. 


T _N^ _ N^ CPP^ 


(5.4) 
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Figure 27. WD for Frank Code with Phase Length =16, CPP = 1, Signal Only, (No.25) 
(a) 2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 28. WD for Frank Code with Phase Length = 16, CPP = 5, Signal Only, (No.28) 
(a) 2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 29. WD for Frank Code with Phase Length = 64, CPP = 1, Signal Only, (No.31) 
(a) 2D Mesh in Frequency Domain, (b) Contour. 
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Figure 30. WD for Fra nk Code with Phase Length = 64, CPP = 5, Signal Only, (No.34) 
(a) 2D Mesh in Frequeney Domain, (b) Contour. 
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2, Summary 


Applying the Wigner Distribution, one ean quickly know the carrier frequency 
(/^) of the Fra nk codes. The carrier frequency is always the first parameter need to be 
identified. Once the carrier frequency is found, by the bandwidth can determine the CPP 
; by the code period or the slope can determine the phase length( ). 


In Figure 28(b) (A^=16, CPP=5) only bandwidth is obtained but not the code 
period. The reason is that the sample number {#sample) is not enough. From (5.3), the 

code period is T = = = 0.080 sec , hence only when > p CQjje 

5 200 fs 


scifnT)l0 512 

period can be shown in the figure. In this chapter, the-—— = — = 0.073 < 0.080. 




7000 


Table 7 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and phase length from Figure 27 to 30. They are done by the 
visible inspection of the WD results. 
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Frank 

No. 

Carrier Frequency 

Bandwidth 

Code Period 

Phase Fength 

25 

100% 

100% 

100% 

100% 

26 

90% 

110% 

100% 

100% 

27 

110% 

110% 

100% 

100% 

28 

95% 

100% 

0% 

0% 

29 

98% 

95% 

0% 

0% 

30 

0% 

0% 

0% 

0% 

31 

100% 

100% 

100% 

100% 

32 

95% 

105% 

100% 

100% 

33 

98% 

105% 

100% 

100% 

34 

100% 

100% 

0% 

0% 

35 

100% 

0% 

0% 

0% 

36 

100% 

0% 

0% 

0% 


Frank code effectiveness on detection 



Parameters 


□ Only signal 
OO dB 

□ (-)6dB 


Table 7. WD Detection Effectiveness for Fra nk Code Signals. 
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c. 


PI CODE 


In case of a double sideband detection (local oscillator is at band center) of a step 
approximation of a linear frequeney modulation, whieh is the PI eode. The PI eode also 
eonsists of elements. If i is the number of the sample in a given frequeney and j is 
the number of the frequeney, the phase of the zth sample ofy'th frequency is [Ref. 20] 

=^[A^-(27-l)][(7-l)iV + (/-l)] (5.5) 

where i = \:N and j = \:N 


Twelve PI eode signals are introdueed in this ehapter. They are numbered from 
37 to 48 and listed in Table 8. A IkHz earner frequency and 7kHz sampling frequency 
are used for these signals. 


No. 

PI 

N (Phase) 

Cyeles/Phase 

SNR 

37 

PI 1 7 4 1 s 

4 

1 

Signal Only 

38 

PI 1 7 4 1 0 

4 

1 

OdB 

39 

PI 1 7 4 1 -6 

4 

1 

-6dB 

40 

PI 1 7 4 5 s 

4 

5 

Signal Only 

41 

PI 1 7 4 5 0 

4 

5 

OdB 

42 

PI 1 7 4 5 -6 

4 

5 

-6dB 

43 

PI 1 7 8 1 s 

8 

1 

Signal Only 

44 

PI 1 7 8 1 0 

8 

1 

OdB 

45 

PI 1 7 8 1 -6 

8 

1 

-6dB 

46 

PI 1 7 8 5 s 

8 

5 

Signal Only 

47 

PI 1 7 8 5 0 

8 

5 

OdB 

48 

PI 1 7 8 5 -6 

8 

5 

-6dB 


Table 8. PI Code Signals. 
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1, Wigner Distribution for PI Code 


Consider the same twelve PI eode signals in Table 8 and use them as the inputs 
for Wigner Distribution (WD). The WD results of these twelve signals are showed in 
Figure 31 to 34. The mesh plots show the frequeney domain of the PI eode signals after 
WD. The contour plots show both the frequency domain and time domain of the results. 


Take Figure 31 as an example. First, the carrier frequency can be clearly 

found by the location of the highest or lowest peak value in Figure 31(a). Secondly, in 
Figure 31(b), the 3dB bandwidth can be measured as lOOOHz. Since 


B=—=^^ 

h cpp 


(5.6) 


so, the cycles per phase (CPP) can be computed by CPP = 


A 

B 


1000 

1000 


For PI code 


7 

T = N%=— (5.7) 

since the code period (T) is measured as 0.016 in Figure 31(b), so the phase length (N^) 
can be calculated from = B-T = 1000• 0.016 = 16 . The slope which measured in 
Figure 31 (b) as 1.6x10^* sec^ can be also use to determine the and CPP by the 
following equation. 

^_T _N^ _N^-CPP^ 
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(b) 

Figure 31. WD for PI Code with Phase Length =16, CPP = 1, Signal Only, (No.37) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 32. WD for PI Code with Phase Length = 16, CPP = 5, Signal Only, (No.40) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 33. WD for PI Code with Phase Length = 64, CPP = 1, Signal Only, (No.43) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 34. WD for PI Code with Phase Length = 64, CPP = 5, Signal Only, (No.46) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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2, Summary 


Applying the Wigner Distribution, one ean quickly know the carrier frequency 
of the Pi codes. The carrier frequency is always the first parameter need to be 
identified. Once the carrier frequency is found, by the bandwidth can determine the CPP 
; by the code period or the slope can determine the phase length( ). 


In Figure 32(b) (A^=16, CPP=5) only bandwidth is obtained but not the code 
period. The reason is that the sample number {#sample) is not enough. From (5.7), the 

code period is T = = = 0.080 sec , hence only when > p CQjje 

5 200 fs 


scifnT)l0 512 

period can be shown in the figure. In this chapter, the-—— = — = 0.073 < 0.080. 




7000 


Table 9 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and phase length from Figure 31 to 34. They are done by the 
visible inspection of the WD results. 
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PI 

No. 

Carrier Frequency 

Bandwidth 

Code Period 

Phase Length 

37 

120% 

70% 

100% 

100% 

38 

110% 

80% 

100% 

100% 

39 

0% 

0% 

0% 

0% 

40 

100% 

100% 

0% 

0% 

41 

100% 

90% 

0% 

0% 

42 

0% 

0% 

0% 

0% 

43 

90% 

50% 

100% 

100% 

44 

90% 

60% 

100% 

100% 

45 

0% 

0% 

0% 

0% 

46 

100% 

100% 

0% 

0% 

47 

100% 

120% 

0% 

0% 

48 

0% 

0% 

0% 

0% 


P1 code effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Table 9. WD Detection Effectiveness for PI Code Signals. 
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D, P2 CODE 


The P2 code also consists elements. If i is the number of the sample in a 
given frequency and j is the number of the frequency, the phase of the zth sample ofyth 
frequency is [Ref 20] 




(5.9) 


where i = l:N and j = 1:N 


Twelve P2 code signals are introduced in this chapter. They are numbered from 


49 to 60 and listed in Table 10. Phase length A^^=16, 64 and CPP=\, 5 are used in this 
chapter. A IkHz carrier frequency and 7kHz sampling frequency are used for these 
signals. 


No. 

P2 

N (Phase) 

Cycles/Phase 

SNR 

49 

P2 1 7 4 1 s 

4 

1 

Signal Only 

50 

P2 1 7 4 1 0 

4 

1 

OdB 

51 

P2 1 7 4 1 -6 

4 

1 

-6dB 

52 

P2 1 7 4 5 s 

4 

5 

Signal Only 

53 

P2 1 7 4 5 0 

4 

5 

OdB 

54 

P2 1 7 4 5 -6 

4 

5 

-6dB 

55 

P2 1 7 8 1 s 

8 

1 

Signal Only 

56 

P2 1 7 8 1 0 

8 

1 

OdB 

57 

P2 1 7 8 1 -6 

8 

1 

-6dB 

58 

P2 1 7 8 5 s 

8 

5 

Signal Only 

59 

P2 1 7 8 5 0 

8 

5 

OdB 

60 

P2 1 7 8 5 -6 

8 

5 

-6dB 


Table 10. P2 Code Signals. 
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1, Wigner Distribution for P2 Code 


Consider the same twelve P2 code signals in Table 10 and use them as the inputs 
for Wigner Distribution (WD). The WD results of these twelve signals are showed in 
Figure 35 to 38. The mesh plots show the frequency domain of the P2 code signals after 
WD. The contour plots show both the frequency domain and time domain of the results. 

Take Figure 35 as an example. First, the carrier frequency can be clearly 

found by the location of the highest or lowest peak value in Figure 35(a). Secondly, in 
Figure 35(b), the 3dB bandwidth can be measured as lOOOHz. Since 

1 f 

B = - = -^ (5.10) 

t, CPP 

so, the cycles per phase {CPP) can be computed by CPP ^ 


T = N\^— (5.11) 

since the code period (7) is measured as 0.016 in Figure 35(b), so the phase length (N^) 
can be calculated from = B-T = 1000• 0.016 = 16 . The slope which measured in 
Figure 35(b) as 1.6x10^* sec^ can be also use to determine the and CPP by the 
following equation. 


T _N^ _N^ CPP^ 


(5.12) 
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Figure 35. WD for P2 Code with Phase Length =16, CPP = 1, Signal Only, (No.49) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 36. WD for P2 Code with Phase Length = 16, CPP = 5, Signal Only, (No.52) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 37. WD for P2 Code with Phase Length = 64, CPP = 1, Signal Only, (No.55) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 38. WD for P2 Code with Phase Length = 64, CPP = 5, Signal Only, (No.58) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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2, Summary 


Applying the Wigner Distribution, one ean quickly know the carrier frequency 
of the P2 codes. The carrier frequency is always the first parameter need to be 
identified. Once the carrier frequency is found, by the bandwidth can determine the CPP 
; by the code period or the slope can determine the phase length( ). 


In Figure 36(b) (A^=16, CPP=5) only bandwidth is obtained but not the code 
period. The reason is that the sample number {#sample) is not enough. From (5.11), the 

code period is T = = = 0.080 sec , hence only when > p CQjje 

5 200 fs 


scifnT)l0 512 

period can be shown in the figure. In this chapter, the-—— = — = 0.073 < 0.080. 




7000 


Table 11 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and phase length from Figure 35 to 38. They are done by the 
visible inspection of the WD results. 
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P2 

No. 

Carrier Frequency 

Bandwidth 

Code Period 

Phase Length 

49 

100% 

70% 

100% 

100% 

50 

110% 

60% 

90% 

100% 

51 

0% 

0% 

0% 

0% 

52 

100% 

100% 

0% 

0% 

53 

100% 

0% 

0% 

0% 

54 

95% 

0% 

0% 

0% 

55 

100% 

100% 

100% 

100% 

56 

100% 

110% 

110% 

100% 

57 

0% 

110% 

110% 

100% 

58 

100% 

100% 

0% 

0% 

59 

100% 

90% 

0% 

0% 

60 

120% 

0% 

0% 

0% 


P2 code effectiveness on detection 


150% 

130% 

110 % 


I 90% 

nj 

O 70% 

o 

Q. 


50% 


30% 


10 % 


□ Only signal 

□ OdB 

□ (-)6dB 


Carrier freq.(Hz) Bandwidth (Hz) Code period (ms) 

Parameters 


Phases 


Table 11. WD Detection Effectiveness for P2 Code Signals. 
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E, 


P3 CODE 


In the case of linear frequency modulation waveforms, the conceptual coherent 
detection using single sideband detection and sampling at the Nyquist rate yield a 
polyphase code call the P3. The phase of the zth sample of the P3 code is given by [Ref. 
20 ] 


where i = \:N 


(5.13) 


Twelve P3 code signals are introduced in this chapter. They are numbered from 
61 to 72 and listed in Table 12. In this chapter, the phase length = 16 and 64 and 
cycles per phase (CPP) = 1 and 5 are used. A IkHz carrier frequency and 7kHz sampling 
frequency are used for these signals. 


No. 

P3 

N (Phase) 

Cycles/Phase 

SNR 

61 

P3 1 7 16 1 s 

16 

1 

Signal Only 

62 

P3 1 7 16 1 0 

16 

1 

OdB 

63 

P3 1 7 16 1 -6 

16 

1 

-6dB 

64 

P3 1 7 16 5 s 

16 

5 

Signal Only 

65 

P3 1 7 16 5 0 

16 

5 

OdB 

66 

P3 1 7 16 5 -6 

16 

5 

-6dB 

67 

P3 1 7 64 1 s 

64 

1 

Signal Only 

68 

P3 1 7 64 1 0 

64 

1 

OdB 

69 

P3 1 7 64 1 -6 

64 

1 

-6dB 

70 

P3 1 7 64 5 s 

64 

5 

Signal Only 

71 

P3 1 7 64 5 0 

64 

5 

OdB 

72 

P3 1 7 64 5 -6 

64 

5 

-6dB 


Table 12. P3 Code Signals. 
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1, Wigner Distribution for P3 Code 


Consider the same twelve P3 code signals in Table 12 and use them as the inputs 
for Wigner Distribution (WD). The WD results of these twelve signals are showed in 
Figure 39 to 42. The mesh plots show the frequency domain of the P3 code signals after 
WD. The contour plots show both the frequency domain and time domain of the results. 

Take Figure 39 as an example. First, the carrier frequency can be clearly 

found by the location of the highest or lowest peak value in Figure 39(a). Secondly, in 
Figure 39(b), the 3dB bandwidth can be measured as lOOOHz. Since 

1 f 

B = - = ^^ (5.14) 

t, CPP 

so, the cycles per phase {CPP) can be computed by CPP ^ 


T = Nt,=j (5.15) 

since the code period (7) is measured as 0.016 in Figure 39(b), so the phase length (N) 
can be calculated from N = B -T = 1000• 0.016 = 16. The slope which measured in 
Figure 39(b) as 1.6x10^* sec^ can be also use to determine the and CPP by the 
following equation. 


T _ N _ N CPP^ 

// 


(5.16) 
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Figure 39. WD for P3 Code with Phase Length = 16, CPP = 1, Signal Only, (No.61) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 40. WD for P3 Code with Phase Length = 16, CPP = 5, Signal Only, (No.64) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 41. WD for P3 Code with Phase Length = 64, CPP = 1, Signal Only, (No.67) (a) 

2D Mesh in Frequency Domain, (b) Contour. 
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Figure 42. WD for P3 Code with Phase Length = 64, CPP = 5, Signal Only, (No.70) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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2, Summary 


Applying the Wigner Distribution, one ean quickly know the carrier frequency 
of the P3 codes. The carrier frequency is always the first parameter need to be 
identified. Once the carrier frequency is found, by the bandwidth can determine the CPP 
; by the code period or the slope can determine the phase length( N). 


In Figure 40(b) (A =16, CPP=5) only bandwidth is obtained but not the code 
period. The reason is that the sample number {#sample) is not enough. From (5.15), the 

code period is T = — = = 0.080sec, hence only when > T , the code period 

5 200 L 


^ safTiDi^ 312 

can be shown in the figure. In this chapter, the-^-=-= 0.073 < 0.080 . 


X 


7000 


Table 13 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and phase length from Figure 39 to 42. They are done by the 
visible inspection of the WD results. 
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P3 

No. 

Carrier Frequency 

Bandwidth 

Code Period 

Phase Length 

61 

95% 

100% 

100% 

100% 

62 

95% 

110% 

100% 

100% 

63 

80% 

110% 

110% 

100% 

64 

100% 

100% 

0% 

0% 

65 

100% 

120% 

0% 

0% 

66 

90% 

0% 

0% 

0% 

67 

110% 

100% 

90% 

100% 

68 

110% 

100% 

90% 

100% 

69 

120% 

110% 

90% 

100% 

70 

100% 

100% 

0% 

0% 

71 

100% 

0% 

0% 

0% 

72 

100% 

0% 

0% 

0% 


P3 code effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Table 13. WD Detection Effectiveness for P3 Code Signals. 
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F. P4 CODE 

The P4 code consists of the discrete phase of the linear chirp waveform taken at 
specific time intervals and exhibits the same range Doppler coupling associated with the 
chirp waveform. The peak sidelobe levels are lower than those of the unweighted chirp 
waveform. The code element of the P4 code are given by [Ref. 20] 

(5.17) 

where i = \:N 


Twelve P4 code signals are introduced in this chapter. They are numbered from 
73 to 84 and listed in Table 14. The phase length or phase elements N are also 16 or 64; 
CPP is 1 or 5 for P4 codes in this chapter. A IkHz carrier frequency and 7kHz sampling 
frequency are used for these signals. 


No. 

P4 

N (Phase) 

Cycles/Phase 

SNR 

73 

P4 1 7 16 1 s 

16 

1 

Signal Only 

74 

P4 1 7 16 1 0 

16 

1 

OdB 

75 

P4 1 7 16 1 -6 

16 

1 

-6dB 

76 

P4 1 7 16 5 s 

16 

5 

Signal Only 

77 

P4 1 7 16 5 0 

16 

5 

OdB 

78 

P4 1 7 16 5 -6 

16 

5 

-6dB 

79 

P4 1 7 64 1 s 

64 

1 

Signal Only 

80 

P4 1 7 64 1 0 

64 

1 

OdB 

81 

P4 1 7 64 1 -6 

64 

1 

-6dB 

82 

P4 1 7 64 5 s 

64 

5 

Signal Only 

83 

P4 1 7 64 5 0 

64 

5 

OdB 

84 

P4 1 7 64 5 -6 

64 

5 

-6dB 


Table 14. P4 Code Signals. 
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1, Wigner Distribution for P4 Code 


Consider the same twelve P4 code signals in Table 14 and use them as the inputs. 
The WD results of these twelve signals are showed in Figure 43 to 46. The mesh plots 
show the frequency domain of the P4 code signals after WD. The contour plots show 
both the frequency domain and time domain of the results. 

Take Figure 43 as an example. First, the carrier frequency can be clearly 

found by the location of the highest or lowest peak value in Figure 43(a). Secondly, in 
Figure 43 (b), the 3dB bandwidth can be measured as lOOOHz. Since 

1 f 

B = - = -^ (5.18) 

t, CPP 

so, the cycles per phase {CPP) can be computed by CPP ^ 


T = Nt, 


B 


(5.19) 


since the code period (T) is measured as 0.016 in Figure 43 (b), so the phase length (A^) 
can be calculated from N = B -T = 1000• 0.016 = 16. The slope which measured in 
Figure 43 (b) as 1.6 x 10sec^ can be also use to determine the and CPP by the 
following equation. 


T _ N _ N CPP^ 

// 


(5.20) 
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(b) 

Figure 43. WD for P4 Code with Phase Length =16, CPP = 1, Signal Only, (No.73) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 44. WD for P4 Code with Phase Length = 16, CPP = 5, Signal Only, (No.76) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 45. WD for P4 Code with Phase Length = 64, CPP = 1, Signal Only, (No.79) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 


90 































































mesh of P 4 ^ 7545^ 


4001- 

g 200 

O) 

T3 

3 

E 0 

O) 

(0 

s 

■200 

-400 

-600 


0.07 - 

0 06 - 

005 

2 0 04 h 
0) 

E 

“ 0 03 
0 02 - 

0.01 - 


500 


■ ■ ..■ 

-■ X!. 





1 : 



















. . i" 

i 





i i 


1 




1000 


1500 2000 

Frequency (Hz) 

(a) 

contour ofP4^7g45^ 


2500 


3000 


Bandwidth (B) = 200 Hz 


I:' 


500 1000 1500 2000 

Frequency (Hz) 


2500 


3000 


(b) 

Figure 46. WD for P4 Code with Phase Length = 64, CPP = 5, Signal Only, (No.82) (a) 

2D Mesh in Frequeney Domain, (b) Contour. 
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2, Summary 


Applying the Wigner Distribution, one ean quickly know the carrier frequency 
of the P4 codes. The carrier frequency is always the first parameter need to be 

identified. Once the carrier frequency is found, by the bandwidth can determine the CPP 
; by the code period or the slope can determine the phase length( N). 


In Figure 44(b) (A =16, CPP=5) only bandwidth is obtained but not the code 
period. The reason is that the sample number {#sample) is not enough. From (5.19), the 

code period is T = — = = 0.080sec, hence only when > T , the code period 

5 200 L 


^ safTiDi^ 312 

can be shown in the figure. In this chapter, the-^-=-= 0.073 < 0.080 . 


X 


7000 


Table 9 below shows the detection effectiveness of the carrier frequency, 3dB 
bandwidth, code period and phase length from Figure 43 to 46. They are done by the 
visible inspection of the WD results. 
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P4 

No. 

Carrier Frequency 

Bandwidth 

Code Period 

Phase Length 

73 

90% 

100% 

100% 

100% 

74 

90% 

95% 

110% 

100% 

75 

80% 

90% 

110% 

100% 

76 

100% 

100% 

0% 

0% 

77 

95% 

0% 

0% 

0% 

78 

100% 

0% 

0% 

0% 

79 

110% 

100% 

100% 

100% 

80 

115% 

105% 

95% 

100% 

81 

125% 

110% 

90% 

100% 

82 

100% 

100% 

0% 

0% 

83 

100% 

110% 

0% 

0% 

84 

100% 

120% 

0% 

0% 


P4 code effectiveness on detection 



□ Only signal 
nOdB 

□ (-)6dB 


Parameters 


Table 15. WD Detection Effectiveness for P4 Code Signals. 
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G. 


COMPARISON OF POLYPHASE CODES 


Table 16 shows the phases for all the different polyphase eodes diseussed in this 

thesis. 


Frank Code 

'2.71 

=-^(z- 1 )( 7 -1) where i = \\N and 7 = 1;A 

PI Code 

A, = (2J - !)][(,• - DA^ + (,• - 1)] .l,ere i = and J = 

P2 Code 

Jwhere z = 1: A and 7 = 1: A 

P3 Code 

61 = — (z - l)(z -1) where z = 1: A 

N 

P4 Code 

( * 1 

6 , = n - nii - 1) where z = 1: A 

A 


Table 16. Phases of All the Polyphase Code. 


Figure 47 compares the WD of all the polyphase codes with same CPP = 1 and 
phase length =16. Figure 48 is the zoom in figure for Figure 47. Figure 49 compares the 
WD of all the polyphase codes with same CPP = 5 and phase length =16. Figure 50 is 
the zoom in figure for Figure 49. Figure 51 compares the WD of all the polyphase codes 
with same CPP = 1 and phase length = 64. Figure 52 is the zoom in figure for Figure 51. 
Figure 53 compares the WD of all the polyphase codes with same CPP = 5 and phase 
length = 64. Figure 54 is the zoom in figure for Figure 53. 
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Figure 47. Wigner Distribution for Polyphase Codes with CPP = 1 and Phase Length ■ 
16 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code, (e) P4 Code. 
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Figure 48. Zoom In for Figure 47 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code. 

(e) P4 Code. 
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Figure 49. Wigner Distribution for Polyphase Codes with CPP = 5 and Phase Length = 
16 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code, (e) P4 Code. 
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Figure 50. Zoom In for Figure 49 (a) Frank Code, (b) PI Code, (e) P2 Code, (d) P3 Code. 

(e) P4 Code. 
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Figure 51. Wigner Distribution for Polyphase Codes with CPP = 1 and Phase Length ■ 
64 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code, (e) P4 Code. 
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Figure 52. Zoom In for Figure 51 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code. 

(e) P4 Code. 
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Figure 53. Wigner Distribution for Polyphase Codes with CPP = 5 and Phase Length ■ 
64 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code, (e) P4 Code. 
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Figure 54. Zoom In for Figure 53 (a) Frank Code, (b) PI Code, (c) P2 Code, (d) P3 Code. 

(e) P4 Code. 


102 










































Some features are observed from the above eight figures. First, for all the 
polyphase signals with CPP =1 (Figures 47, 48, 51 and 52), the P2 codes have negative 
slope, which is different from other polyphase signals. Secondly, for the WD analysis, 
when CPP =1 (Figures 47, 48, 51 and 52), the PI code is the most difficult to determine 
the 3 dB bandwidth as well as the code period. On the contrary, the P3 and P4 code can 
provide very clear parameters in Figure 48 and 52. Thirdly, when comparing with all the 
zoom in figures (Figures 48, 50, 52 and 54), the P3 code and the P4 code have the best 
parameter resolution among the five polyphase signals. 
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VI. COSTAS FREQUENCY HOPPING 


A, COSTAS 

A Costas array is an nxn array of frequencies and times with exactly one 
frequency in each row and column, and with distinct vector differences between all pairs 
of frequencies. As a frequency-hopping pattern for radar or sonar, a Costas array has an 
optimum ambiguity function, since any translation of the array parallel to the coordinate 
axes produces at most one out of phase coincidence. [Ref 21] 

In a frequency hopping system, the signal consists of one or more frequencies 
being chosen from a set available frequencies, for transmission at each 

of a set of consecutive time intervals. For modeling purposes, it is reasonable 

to consider the situation in which m = n, and a different one of n equally spaced 
frequencies {yi,/ 2 ,••••/„} transmitted during each of the n equal duration time intervals 

• Such a signal is represented by a nxn permutation matrix A, where the n 
rows correspond to the n frequencies, the n columns correspond to the n intervals, and the 
entry a^j equals 1 means transmission and 0 otherwise. 

For example if the coding sequence {a-} is {4,7,1,6,5,2,3,}, then the coding 

matrix, the difference matrix and the ambiguity sidelobe matrix of a Costas signal are 
shown as follows. [Ref. 22] 
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Figure 55. The Coding Matrix, Difference Matrix and Ambiguity Sidelobe Matrix of a 

Costas Signal. 


The element of the difference matrix in row i and column j, is 


D,-=a-^j-aj, i + j<N (6.1) 

with blanks at the remaining locations. Thus the first row is formed by taking differences 
between adjacent terms in the sequence, the second row by taking differences between 
next adjacent terms, and so on. The corresponding values in the difference matrix and 
sidelobe matrix were marker with different color circles. [Ref. 22] 


The difference matrix not only determines whether a sequence is Costas or not, 
but it also identifies the major sidelobe locations in the positive delay slots of the 
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ambiguity function. The sidelobes in the negative delay slots are obtained by applying 


the rule of symmetry with respeet with the origin. [Ref 22] 


Twelve Costas signals are introdueed in this ehapter. They are numbered from 85 
to 96 and listed in Table 17. The “Sequenee” means the frequency hopping sequence. 
There are two different sequences used in this table. They are the carrier frequeney in 
order of the number in the sequenee. The unit of these carrier frequeneies is kHz. The 
“Cycles/ Frequency” means the number of eycles per frequency, the number 10 and 20 
are used here. The sampling frequency for these Costas signals is 15kHz for sequence 
4716523 and 17kHz for sequence 2638751. 


No. 

COSTAS 

Sequence 

Cy e les/F requeney 

SNR 

85 

C 1 15 10 s 

4716523 

10 

Signal Only 

86 

C 1 15 10 0 

4716523 

10 

OdB 

87 

C 1 15 10 -6 

4716523 

10 

-6dB 

88 

C 1 15 20 s 

4716523 

20 

Signal Only 

89 

C 1 15 20 0 

4716523 

20 

OdB 

90 

C 1 15 20 -6 

4716523 

20 

-6dB 

91 

C 2 17 10 s 

2638751 

10 

Signal Only 

92 

C 2 17 10 0 

2638751 

10 

OdB 

93 

C 2 17 10 -6 

2638751 

10 

-6dB 

94 

C 2 17 20 s 

2638751 

20 

Signal Only 

95 

C 2 17 20 0 

2638751 

20 

OdB 

96 

C 2 17 20 -6 

2638751 

20 

-6dB 


Table 17. COSTAS Signals. 
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B, 


WIGNER DISTRIBUTION FOR COSTAS 


Consider the same twelve Costas frequeney-hopping signals in Table 17 and use 
them as the inputs for Wigner Distribution (WD). The WD results of signals number 85, 
88, 91, 94 are showed in Figures 56 to 59. The mesh plots show the frequeney domain of 
the Costas frequency-hopping signals after WD. The contour plots show both the 
frequency domain and time domain of the results. 

First, the seven carrier frequencies and the entire cross terms are shown in every 
figure. Since the carrier frequencies are IkHz, 2kHz, 3kHz, 4kHz, 5kHz, 6kHz and 
7kHz, the cross terms will be 

|{(l+2),(l+3),(l+4),(l+5),(l+6),(l+7) 

,(2+3),(2+4),(2+5),(2+6),(2+7) 

,(3+4),(3+5),(3+6),(3+7) 

,(4+5),(4+6),(4+7) 

,(5+6),(5+7) 

,(6+7)} kHz 

From the above calculation, the cross terms appear at 1.5 kHz one time, 2 kHz one time, 
2.5 kHz two times, 3 kHz two times, 3.5 kHz three times, 4 kHz three times, 4.5 kHz 
three times, 5 kHz two times, 5.5 kHz two times, 6 kHz one time, 6.5 kHz one time. In 
Figure 349 to 360, only the IkHz and 7kHz are real carrier frequencies, all the others are 
possibly carrier frequencies or cross terms. 

Secondly, in Figure 56(b) to 59(b) there marked with the blue arrows, which 

indicate the frequency sequence from the top to the bottom. Take Figure 56(b) as an 

example, the carrier frequency sequence is, as the blue arrows shown from the top to the 

bottom, AkHz IkHz IkHz 6kHz 5kHz 2kHz 3kHz. 
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As can be seen, the eolors of earner frequencies and eross terms in Figure 57(b) 
are different. The eolor of the cross terms is red line enclosed by yellow, but the real 
carrier frequeneies display different eolor and shapes. They are yellow cireles with green 
curves on the right and left hand sides. Aceording to Figure 10, the frequeney of the 
eross terms is half of the two real earner frequencies. The time period of the eross terms 
is also half of the time period sum of the two real carrier frequencies. 

Taking a look of the 4kHz in Figure 57(b), there are four 4kHz bars on this plot. 
From the top to the bottom, the first yellow bar is the real carrier frequeney. The second 
red bar is a cross term eaused by 7kHz and IkHz, the appearing time period is right on 
the eenter of those of in 7kHz and IkHz. The third red bar is a eross term eaused by 
6kHz and 2kHz, and with time period in the center of those in 6kHz and 2kHz. The forth 
red bar is a eross term caused by 5kHz and 3kHz with time period in the eenter of those 
in 5kHz and 3kHz. Note there are two 3kHz red bars in Figure 57(b). The first is caused 
by 4kHz and 2kHz carrier frequencies, the second one is caused by IkHz and 5kHz 
carrier frequeneies. Another 3kHz yellow bar is a real carrier frequeney. Figure 58 and 
Figure 59 show another sequence, the earrier frequeneies are also examined in the 
figures. 


In these figures the time scale depends on the sampling frequency and the number 


sofyiDip 

of sample points. The total time is-, whieh can be used to determine the time in 

f s 


eaeh frequency. Since the total sample number is too big to use the WD, in this ehapter 
WD proeesses 350 points for Cyeles/Frequeney = 10 and 700 points for 
Cyeles/Frequeney = 20. 
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Figure 56. WD for COSTAS Frequency Flopping, Signal Only, (No.85) (a) 2D Mesh in 

Frequency Domain, (b) Contour. 
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Figure 57. WD for COSTAS Frequency Flopping, Signal Only, (No.88) (a) 2D Mesh in 

Frequency Domain, (b) Contour. 
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Figure 58. WD for COSTAS Frequency Flopping, Signal Only, (No.91) (a) 2D Mesh in 

Frequency Domain, (b) Contour. 
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Figure 59. WD for COSTAS Frequency Flopping, Signal Only, (No.94) (a) 2D Mesh in 

Frequency Domain, (b) Contour. 
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c. 


SUMMARY 


To analyze a signal using Wigner distribution, the eross terms are unavoidable for 
multi-frequencies as well as frequency-hopping signals, like Costas. The cross terms 
have the same magnitude as the real carrier frequencies. The frequency of any cross term 
is half of the sum of any two real carrier frequencies. For example, if and f ^2 


carrier frequencies, then the cross term generated by these two frequencies is 



Both real carrier frequencies and cross term exist in the mesh plots and contour 
plots. Fortunately, the contour plots still can see the sequence of the Costas frequency 
hopping signals. As long as the carrier frequencies and their sequence have been found, 
then the COSTAS frequency hopping signals can be recovered. 

Table 18 below shows the detection effectiveness of the carrier frequency 
sequence, the time in each sequence and the code period from Figure 56 to 59. They are 
done by the visible inspection of the WD results. 
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COSTAS 

No. 

Sequence 

Time in Each Erequency 

Code Period 

85 

100% 

100% 

100% 

86 

80% 

50% 

50% 

87 

0% 

0% 

0% 

88 

100% 

100% 

100% 

89 

80% 

70% 

70% 

90 

0% 

0% 

0% 

91 

100% 

100% 

100% 

92 

80% 

50% 

50% 

93 

0% 

0% 

0% 

94 

100% 

100% 

100% 

95 

80% 

70% 

70% 

96 

0% 

0% 

0% 


Costas code effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 

□ (-)6dB 


Table 18. WD Detection Effectiveness for COSTAS Code Signals. 
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VII. PHASE SHIFT KEYING/FREQUENCY SHIFT KEYING 

(PSK/FSK) 

A, PSK/FSK USING A COSTAS-BASED FREQUENCY-HOPPING 

TECHNIQUE 

The PSK/FSK using Costas modulation technique is the result of a combination of 
frequency-shift keying based on a Costas frequency-hopping matrix and phase-shift 
keying using Barker sequences with different lengths. Here will briefly describe the 
phase encoding applied to a Costas signal, generating the FSK/PSK combined waveform. 

Consider a Costas frequency-hopped signal, the firing order defines which 
frequencies will appear and with which duration. For CW radars, the usual terminology 
does not apply to this case. Instead of a “bursf’ of pulses, frequencies have being 
continuously emitted during a defined period of time. This period may be divided into 
sub-periods, labeled Tp for each frequency. The length of each sub-period depends on 
the sampling interval. During each sub-period, as the signal stays in one of the 
frequencies, a binary phase modulation occurs according to a Barker sequence of length 

five (+ -h -h - -h), seven (+ -i- -i- - - -i- -), eleven (+ + +- -i----i--)or thirteen 

+ + - + - +). [Ref. 23] 

For example, the FSK/PSK signal defined by: S = {V, V, V, F, V, 2"^, 2^, 2"^, 2', 
2^, 3^, 3"^, 3^, 3', 3^, 4"^, 4^, 4^, 4', 4^, 5"^, 5^, 5^, 5', 5"^} represents a waveform comprised 
of five different frequencies, which are each subdivided into five phase slots, labeled Tp , 
according to the Barker sequence of length five (+ -I- -I- - -h). 
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The final waveform may be seen as a binary phase-shifting modulation within 
eaeh frequency hop, resulting in 25 phase slots equally distributed in each frequency. 

If consider as the number of frequency hops and Np as the number of phase 
slots of duration Tp in each frequency sub-period Tp , the total number of phase slots in 
the FSK/PSK waveform is given by: N = Np ■ Np. [Ref 23] 


The Barker sequence is generated and the frequency-hopping signal is then phase- 
modulated accordingly. For example, if the first Costas sequence is selected, after a 
phase modulation using a Barker sequence of length 5, the final waveform becomes: 


S = 4^, 4^, 4^, 4', 4^, 7^, 7^, 7^, 7', 7^, 1^, 1^, 1^, 1', 1^, 6^, 6^, 6^, 6', 6^, 5^, 5^, 5^, 


5', 5^ 2+ 2+ 2+ 2', 2^ 3^ 3^ 3^ 3', 3^ [Ref 23] 


B, FSK/ PSK COMBINED USING A TARGET-MATCHED FREQUENCY 
HOPPING 


Instead of spreading the energy of the signal equally over a broad bandwidth, this 
type of technique concentrates the signal energy in specific spectral locations of most 
importance for the radar and its typical targets, within the broad-spectrum bandwidth. 
The produced signals have a pulse compression characteristic, and therefore they can 
achieve a low probability of intercept. 

The implementation starts with a simulated-target time-radar response. This data 
is the Fourier transformed, the frequency components, their correspondent magnitudes, 
and their initial phases. A random selection process chooses each frequency with a 
probability distribution function defined by the spectral characteristics of the target of 

interest obtained from the Fast Fourier Transform so that frequencies at the spectral peaks 
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of the target (highest magnitudes) are transmitted more often. Eaeh “frequency hop,” 
transmitted during a specific period of time, is also modulated in phase, having its initial 
phase value modified by a pseudo-random spreading-phase sequence code of values 
equally likely to be zero or n radians. 

The matched FSK/PSK radar will then use a correlation receiver with a phase 
mismatched reference signal instead of a perfectly phase matched reference. This allows 
the radar to generate signals that can match a target’s spectral response in both magnitude 
and phase. 

This simulation serves the purpose of testing the performance of a new digital 
cyclostationary receiver against these kinds of signals. A better implementation would 
require an iterative solution since the response of a target is not always known 
beforehand. [Ref. 23] 

C. PSK/FSK AND TEST SIGNALS PAF AND PSD ANALYSIS 

Eight PSK/ESK COSTAS, eight PSK/FSK target and two single carrier frequency 
signals are introduced in this chapter. They are numbered from 97 to 114 and listed in 
Table 19. The “Sequence” means the frequency hopping sequence. The numbers in the 
sequence are the order of the carrier frequencies. The unit of these carrier frequencies is 
kHz. The “Cycles/ Phase” means the number of cycles per phase, the number 1, 5 and 10 
are used here. The sampling frequency is 15kHz for these PSK/FSK signals, 7kHz for 
the two test signals. 
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No. 

FSK/PSK COSTAS 

Sequence 

Barker 

Code 

Length 

Cycles/ 

Phase 

SNR 

97 

FSK PSK C 1 15 5 5 s 

4716523 

5 

5 

Signal Only 

98 

FSK PSK C 1 15 5 5 0 

4716523 

5 

5 

OdB 

99 

FSK PSK C 1 15 5 1 s 

4716523 

5 

1 

Signal Only 

100 

FSK PSK C 1 15 5 1 0 

4716523 

5 

1 

OdB 

101 

FSK PSK C 1 15 11 5 s 

4716523 

11 

5 

Signal Only 

102 

FSK PSK C 1 15 11 5 0 

4716523 

11 

5 

OdB 

103 

FSK PSK C 1 15 11 1 s 

4716523 

11 

1 

Signal Only 

104 

FSK PSK C 1 15 11 1 0 

4716523 

11 

1 

OdB 


No. 

FSK/PSK Target 

Sequence 

Random 

Hops 

Cycles/ 

Phase 

SNR 

105 

FSK PSK T 15 128 5 s 

1 

128 

5 

Signal Only 

106 

FSK PSK T 15 128 5 0 

1 

128 

5 

OdB 

107 

FSK PSK T 15 256 5 s 

1 

256 

5 

Signal Only 

108 

FSK PSK T 15 256 5 0 

1 

256 

5 

OdB 

109 

FSK PSK T 15 128 10 s 

2 

128 

10 

Signal Only 

no 

FSK PSK T 15 128 10 0 

2 

128 

10 

OdB 

111 

FSK PSK T 15 256 10 s 

2 

256 

10 

Signal Only 

112 

FSK PSK T 15 256 10 0 

2 

256 

10 

OdB 


No. 

Test Signal 

fc (kHz) 

/, (kHz) 

SNR 

113 

T 1 7 1 s 

1 

7 

Signal Only 

114 

T 2 7 1 s 

1 and 2 

7 

Signal Only 


Table 19. FSK/PSK and Test Signals. 


D, WIGNER DISTRIBUTION FOR PSK/FSK AND TEST SIGNALS 

Consider the same eight PSK/FSK COSTAS, eight PSK/FSK target and two 
single earner frequency signals introduced in Table 19 and use them as the inputs for 
Wigner Distribution (WD). The WD results of these twelve signals are showed in Figure 
60 to 67. The mesh plots show the frequency domain of these signals after WD. The 
contour plots show both the frequency domain and time domain of the results. 

The following figures for the PSK/FSK COSTAS signals are very similar to the 


figures of COSTAS signals. Comparing with Figure 60(b) and Figure 56(b) and 57(b), 
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the Wigner distribution analysis gets very similar results for the PSK/FSK COSTAS 
signals and COSTAS signals. The earrier frequeneies sequenee, 4716523, ean be seen as 
yellow eireles and the frequeney cross terms can be seen as red bars in Figure 60(b). As 
a matter of fact, no matter whether the 5 bits or 11 bits Barker code length, the time- 
frequency relationships are the same for the PSK/FSK COSTAS signals. 

In Figure 61(b) and 63(b), the CPP is one, so the shape of the carrier frequencies 
are different from those with five cycles per phase in Figure 60(b) and 62(b). One can 
still distinguish the real carrier frequencies by the lighter colors. 

Since the phase changes are by random, for the PSK/FSK target signals neither 
the carrier frequencies nor the cross terms can be examined on the mesh plots and 
contour plots. Figure 68 shows the one carrier frequency test signal. Figure 69 shows the 
two carrier frequencies and frequency cross term as described in Chapter II. 
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Figure 60. WD for FSK/PSK COSTAS with 5 bit Barker Code, CPP = 5, Signal Only, 
(No.97) (a) 2D Mesh in Frequency Domain, (b) Contour. 
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Figure 61. WD for FSK/PSK COSTAS with 5 bit Barker Code, CPP = 1, Signal Only, 
(No.99) (a) 2D Mesh in Frequeney Domain, (b) Contour. 
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(b) 

Figure 62. WD for FSK/PSK COSTAS with 11 bit Barker Code, CPP = 5, Signal Only, 
(No. 101) (a) 2D Mesh in Frequency Domain, (b) Contour. 
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(b) 

Figure 63. WD for FSK/PSK COSTAS with 11 bit Barker Code, CPP = 1, Signal Only, 
(No. 103) (a) 2D Mesh in Frequency Domain, (b) Contour. 
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Figure 64. WD for FSK/PSK Target with 128 Random Hops, CPP = 5, Signal Only, 
(No. 105) (a) 2D Mesh in Frequency Domain, (b) Contour. 
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Figure 65. WD for FSK/PSK Target with 256 Random Hops, CPP = 5, Signal Only, 
(No. 107) (a) 2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 66. WD for FSK/PSK Target with 128 Random Hops, CPP = 10, Signal Only, 
(No. 109) (a) 2D Mesh in Frequeney Domain, (b) Contour. 
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Figure 67. WD for FSK/PSK Target with 256 Random Hops, CPP = 10, Signal Only, 
(No.l 11) (a) 2D Mesh in Frequency Domain, (b) Contour. 
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Figure 68. WD for 1 kHz Carrier Frequency Test Signal, (No.l 13) (a) 2D Mesh in 

Frequency Domain, (b) Contour. 
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Mesh of two carrier frequencies example 



(a) 


Contour of two carrier frequencies example 



Figure 69. WD for IkHz and 2 kHz Carrier Frequency Test Signal, (No.l 14) (a) 2D 

Mesh in Frequency Domain, (b) Contour. 
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E, SUMMARY 

From the above Figures, only the PSK/FSK COSTAS signals and the last two test 
signals show the carrier frequencies or their sequence. Wigner Distribution cannot 
analyze the PSK/FSK target frequency hopping signals. 

Table 20 below shows the detection effectiveness of the sequence, bandwidth, 
code period and phase from Figure 60 to 67. 


FSK/PSK_COSTAS 

No. 

Sequence 

Bandwidth 

Code Period 

97 

100% 

100% 

100% 

98 

50% 

50% 

50% 

99 

100% 

100% 

100% 

100 

50% 

50% 

50% 

101 

100% 

100% 

100% 

102 

50% 

50% 

50% 

103 

100% 

100% 

100% 

104 

50% 

50% 

50% 

FSK/PSK_TARGET 

No. 

Sequence 

Bandwidth 

Phase 

105 

0% 

0% 

0% 

106 

0% 

0% 

0% 

107 

0% 

0% 

0% 

108 

0% 

0% 

0% 

109 

0% 

0% 

0% 

no 

0% 

0% 

0% 

111 

0% 

0% 

0% 

112 

0% 

0% 

0% 
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FSK/PSK Costas effectiveness on detection 



□ Only signal 
nOdB 


Parameters 


FSK/PSK Target effectiveness on detection 



Parameters 


□ Only signal 

□ OdB 


Table 20. WD Detection Effectiveness for FSK/PSK Signals. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A, CONCLUSIONS 

Nine types of Low Probability of Intereept (LPI) signals were eonsidered in this 
thesis. They are Binary Phase Shift Keying (BPSK), Frequency Modulation Continuous 
Wave (FMCW), Frank code, PI code, P2 code, P3 code, P4 code, COSTAS frequency 
hopping and Phase Shift Keying/Frequency Shift Keying (PSK/FSK) signals. There are 
114 signals were generated for Wigner distribution analysis in this thesis. The detail 
parameters of these 114 test signals are listed in Appendix B. 

The Wigner distribution is a two-dimension function describing the frequency 
content of a signal as a function of time. The Wigner distribution has been noted as one 
of the more useful time-frequency analysis techniques for signal processing. Using the 
Wigner distribution one can identify the frequency and time changes in most of the LPI 
signals. This thesis successfully examined the nine types of LPI signals mentioned 
above. All the signal parameters excluding the PSK/FSK signals can be represented by 
Wigner distribution. For those detectable signals, some of them can show the parameters 
at both SNR = OdB and -6 dB, some of them cannot be detected at SNR = -6dB. That is 
to say the signal to noise ratio is a main concern to the analysis of the LPI signals. Please 
see the technical report. [Ref. 7] 

Other than the signal to noise ratio, there are some other important characteristics 
that need to be emphasized. For BPSK signals, the number of Barker bits can be 
displayed clearly only when the number of cycles per barker bit is one (NPBB = 1). For 
FMCW signals, the parameters can be displayed meaningfully when the modulation 
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period less than or equal to half of the sample points divided by sampling frequeney 




# sample 
2 /. 


) and when the modulation bandwidth less than or equal to 2 times of 


carrier frequency (AF<2/^) or ( AF <2(fJ2-fj ). For the poly phase signals 
including Frank codes, PI codes, P2 codes, P3 codes and P4 codes, the phase length {N or 
N^) can be obtained precisely when the number of cycles per phase is one (CPP = 1). 
The higher the CPP, the harder it is to get the right phase length. 


Last, the frequency cross terms are irremovable for COSTAS frequency hopping 
signals or any other multi-carrier-frequency signal. The cross terms have the same 
amplitude as real carrier frequencies and their frequencies are the average of any two real 

carrier frequencies 

B. RECOMMENDATIONS 


The Choi-Williams Distribution (CWD) is a time frequency analysis technique 
that evolved from Wigner Distribution (WD). The most important advantage of using the 
Choi-Williams distribution is that it has the ability to eliminate the frequency cross terms. 
The definition equation of the Choi-Williams Distribution was introduced in chapter I, 
but no experiments were examined in this thesis. The CWD is deserving of further study. 

The output of the Wigner distribution is a real two-dimensional matrix with the 
size of length (7) times length (0, where / and Q are the in-phase and the quadrature 
component of the input signal. The input signal in this thesis was defined as x = / + j*Q. 
If the length of / or g is 10k points, then the size of the output matrix will be lOOM 
points. In other words, the output matrix will occupy 800MB memory. 
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Thus, the bigger the input signal size, the better and more information will be 
shown in the mesh and contour plots. But the size of the output matrix is too great for 
MATLAB® processing on the computer used in this study; MATLAB® ran out of 

memory. 

The specifications of the computer used to run the Wigner distribution 
MATLAB® program is Intel Pentium© 4 CPU with 2.0 GHz clock speed and 1 GB 
RAM. It is representative of the best in the personal computer market. Future Wigner 
distribution analysis obviously requires more capable computers. 
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APPENDIX A. MATLAB CODE 


A, WIGNERM.M 


function [W] = wigner_fft(x) 

% Design a variable W as the Wigner distribution function for 
input signal x 

% Paul Milne 

% EC4700 Project in Naval Postgraduate School 
% May 2001 
"0 

% Modified by Jen-Yu, Gau 
% Jan. 2002 

N = length(x)/2; 
input signal x 
odd = 0; 
even 

if mod(2*N,2) ~= 0 % If the signal has odd length 

disp('WARNING: Input Sequence Odd:'); 

disp (' A zero will be appended for calculations.'); 

% Display warning messages 

X = [0 x]; % Add a zero at the beginning to make even length 

N = length(x)/2; % Re-compute length 2N 

odd =1; % Mark the original input signal has odd length 

end 

M = 2*N-1 
a = -N; 
b = N-1; 
y = x; 
signal 

for n = a:b % Set index n from -N to N-1 

X = shiftz(y,n); % Call function shiftz.m to shift the 

initial input signal x 

xl = [x(l,N: M) 0 x(l,l:N-l)]; 

x2 = [fliplr(X(1,1:N)) 0 fliplr(x(1,N+1: M))]; 

% Compute the shifted signal xl and x2 in order to get the 
kernel 

kern = xl.*conj (x2); % Compute the kernel as in Eq(2.5) 

W(n-a+l,:) = fft(kern); 

% Apply the fast Fourier transform (FFT) algorithm of Matlab 
command 

end 

W=2*real(W) %Wigner distribution results are always real 

% The final result of the Wigner distribution as in Eq(2.13) 

if odd % If the original input signal has odd length 


% Set M = 2N-1 

% Set a as integration down limit 
% Set b as integration up limit 

% Set y as another name for the input 


% Set 2N as the length of the 
% Initial the signal length as 
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% Then remove extra output created by the adding zero 
W = W(2:2*N,1:2*N-1); % The final result becomes 

end 


B, SHIFTZ,M 


function [y] = shift(x,n) 

% The function is used to shift input signal x with n zeros to the 
right and left 

% in order to get the kernel as in Eg(2.5) 

% When X length 2N, then n is from -N to N-1 

% Paul Milne 

% EC4700 Naval Postgraduate School 
% May 2001 
"0 

% Modified by Jen-Yu, Gau 
% Jan. 2002 

N = length(x)/2; % Set 2N as the length of the input signal x 

if n >= 0 % First shift the up half part of x 

y = [zeros(l,n) x(l:2*N-n)]; % Shift x to the right with n 

zeros 

else % Then shift the down half part of x 

n = abs(n); 

y = [x(l+n:2*N) zeros(l,n)]; 

% Shift X to the left with n zeros 
end 
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APPENDIX B. TEST MATRIX 


No. 

BPSK 

Barker Code 
Length 

NPBB 

SNR 

1 

B 1 7 7 1 s 

7 

1 

Signal Only 

2 

B 1 7 7 1 0 

7 

1 

OdB 

3 

B 1 7 7 1 -6 

7 

1 

-6dB 

4 

B 1 7 11 1 s 

11 

1 

Signal Only 

5 

B 1 7 11 1 0 

11 

1 

OdB 

6 

B 1 7 11 1-6 

11 

1 

-6dB 

7 

B 1 7 7 5 s 

7 

5 

Signal Only 

8 

B 1 7 7 5 0 

7 

5 

OdB 

9 

B 1 7 7 5 -6 

7 

5 

-6dB 

10 

B 1 7 11 5 s 

11 

5 

Signal Only 

11 

B 1 7 11 5 0 

11 

5 

OdB 

12 

B 1 7 11 5 -6 

11 

5 

-6dB 


FMCW 

Mod. BW 
(Hz) 

Mod. Period 
(ms) 

SNR 

13 

F 1 7 250 20 s 

250 

20 

Signal Only 

14 

F 1 7 250 20 0 

250 

20 

OdB 

15 

F 1 7 250 20 -6 

250 

20 

-6dB 

16 

F 1 7 250 30 s 

250 

30 

Signal Only 

17 

F 1 7 250 30 0 

250 

30 

OdB 

18 

F 1 7 250 30 -6 

250 

30 

-6dB 

19 

F 1 7 500 20 s 

500 

20 

Signal Only 

20 

F 1 7 500 20 0 

500 

20 

OdB 

21 

F 1 7 500 20 -6 

500 

20 

-6dB 

22 

F 1 7 500 30 s 

500 

30 

Signal Only 

23 

F 1 7 500 30 0 

500 

30 

OdB 

24 

F 1 7 500 30 -6 

500 

30 

-6dB 


FRANK 

N (Phase) 

Cycles/Phase 

SNR 

25 

FR 1 7 4 1 s 

4 

1 

Signal Only 

26 

FR 1 7 4 1 0 

4 

1 

OdB 

27 

FR 1 7 4 1 -6 

4 

1 

-6dB 

28 

FR 1 7 4 5 s 

4 

5 

Signal Only 

29 

FR 1 7 4 5 0 

4 

5 

OdB 

30 

FR 1 7 4 5 -6 

4 

5 

-6dB 

31 

FR 1 7 8 1 s 

8 

1 

Signal Only 

32 

FR 1 7 8 1 0 

8 

1 

OdB 

33 

FR 1 7 8 1 -6 

8 

1 

-6dB 

34 

FR 1 7 8 5 s 

8 

5 

Signal Only 

35 

FR 1 7 8 5 0 

8 

5 

OdB 

36 

FR 1 7 8 5 -6 

8 

5 

-6dB 


PI 

N (Phase) 

Cycles/Phase 

SNR 
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37 PI 1 7 4 1 s 

38 PI 1 7 4 1 0 

39 PI 1 7 4 1 -6 

40 PI 1 7 4 5 s 

41 PI 1 7 4 5 0 

42 PI 1 7 4 5 -6 

43 PI 1 7 8 1 s 

44 PI 1 7 8 1 0 

45 PI 1 7 8 1 -6 

46 PI 1 7 8 5 s 

47 PI 1 7 8 5 0 

48 PI 1 7 8 5 -6 

_ P2 

49 P2 1 7 4 1 s 

50 P2 1 7 4 1 0 

51 P2 1 7 4 1 -6 

52 P2 1 7 4 5 s 

53 P2 1 7 4 5 0 

54 P2 1 7 4 5 -6 

55 P2 1 7 8 1 s 

56 P2 1 7 8 1 0 

57 P2 1 7 8 1 -6 

58 P2 1 7 8 5 s 

59 P2 1 7 8 5 0 

60 P2 1 7 8 5 -6 

_ P3 

61 P3 1 7 16 1 s 

62 P3 1 7 16 1 0 

63 P3 1 7 16 1 -6 

64 P3 1 7 16 5 s 

65 P3 1 7 16 5 0 

66 P3 1 7 16 5 -6 

67 P3 1 7 64 1 s 

68 P3 1 7 64 1 0 

69 P3 1 7 64 1 -6 

70 P3 1 7 64 5 s 

71 P3 1 7 64 5 0 

72 P3 1 7 64 5 -6 

_ P4 

73 P4 1 7 16 1 s 

74 P4 1 7 16 1 0 

75 P4 1 7 16 1 -6 

76 P4 1 7 16 5 s 

77 P4 1 7 16 5 0 

78 P4 1 7 16 5 -6 


N (Phase) 

4 

4 

4 

4 

4 

4 

8 

8 

8 

8 

8 

8 

N (Phase) 

16 

16 

16 

16 

16 

16 

64 

64 

64 

64 

64 

64 

N (Phase) 

16 

16 

16 

16 

16 

16 


Cycles/Phase 

1 

1 

1 

5 

5 

5 

1 

1 

1 

5 

5 

5 

Cycles/Phase 

1 

1 

1 

5 

5 

5 

1 

1 

1 

5 

5 

5 

Cycles/Phase 

1 

1 

1 

5 

5 

5 


Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 
SNR 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 
SNR 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 
SNR 

Signal Only 
OdB 
-6dB 

Signal Only 
OdB 
-6dB 













































79 

P4 1 7 64 1 s 

64 

1 

Signal Only 

80 

P4 1 7 64 1 0 

64 

1 

OdB 

81 

P4 1 7 64 1 -6 

64 

1 

-6dB 

82 

P4 1 7 64 5 s 

64 

5 

Signal Only 

83 

P4 1 7 64 5 0 

64 

5 

OdB 

84 

P4 1 7 64 5 -6 

64 

5 

-6dB 


COSTAS 

Sequence 

Cycles/Freq, 

SNR 

85 

C 1 15 10 s 

4716523 

10 

Signal Only 

86 

C 1 15 10 0 

4716523 

10 

OdB 

87 

C 1 15 10 -6 

4716523 

10 

-6dB 

88 

C 1 15 20 s 

4716523 

20 

Signal Only 

89 

C 1 15 20 0 

4716523 

20 

OdB 

90 

C 1 15 20 -6 

4716523 

20 

-6dB 

91 

C 2 17 10 s 

2638751 

10 

Signal Only 

92 

C 2 17 10 0 

2638751 

10 

OdB 

93 

C 2 17 10 -6 

2638751 

10 

-6dB 

94 

C 2 17 20 s 

2638751 

20 

Signal Only 

95 

C 2 17 20 0 

2638751 

20 

OdB 

96 

C 2 17 20 -6 

2638751 

20 

-6dB 


FSK/PSK COSTAS 

Sequence 

Barker 

Code 

Length 

Cycles/ 

Phase 

SNR 

97 

FSK PSK C 1 15 5 5 s 

4716523 

5 

5 

Signal Only 

98 

FSK PSK C 1 15 5 5 0 

4716523 

5 

5 

OdB 

99 

FSK PSK C 1 15 5 1 s 

4716523 

5 

1 

Signal Only 

100 

FSK PSK C 1 15 5 1 0 

4716523 

5 

1 

OdB 

101 

FSK PSK C 1 15 11 5 s 

4716523 

11 

5 

Signal Only 

102 

FSK PSK C 1 15 11 5 0 

4716523 

11 

5 

OdB 

103 

FSK PSK C 1 15 11 1 s 

4716523 

11 

1 

Signal Only 

104 

FSK PSK C 1 15 11 1 0 

4716523 

11 

1 

OdB 


FSK/PSK Target 

Sequence 

Random 

Hops 

Cycles/ 

Phase 

SNR 

105 

FSK PSK T 15 128 5 s 

1 

128 

5 

Signal Only 

106 

FSK PSK T 15 128 5 0 

1 

128 

5 

OdB 

107 

FSK PSK T 15 256 5 s 

1 

256 

5 

Signal Only 

108 

FSK PSK T 15 256 5 0 

1 

256 

5 

OdB 

109 

FSK PSK T 15 128 10 s 

2 

128 

10 

Signal Only 

no 

FSK PSK T 15 128 10 0 

2 

128 

10 

OdB 

111 

FSK PSK T 15 256 10 s 

2 

256 

10 

Signal Only 

112 

FSK PSK T 15 256 10 0 

2 

256 

10 

OdB 


Test Signal 

/c(kHz) 

/.(kHz) 

SNR 

113 

T 1 7 1 s 

1 

7 

Signal Only 

114 

T 2 7 1 s 

1 and 2 

7 

Signal Only 
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